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ABSTRACT 


the  nucleate  pool-boiling  performance 
il  mixtures  from  enhanced  evaporator 
tube  surfaces,  an  experimental  apparatus  was  designed, 
constructed  and  instrumented.  The  evaporator  was  made  of  a 
T-shaped  Pyrex  glass  container.  Boiling  occurred  from  a 
smooth,  hard-copper  tube,  15.9  mm  in. )  in  outer  diam¬ 

eter,  12.7  mm  (1/2  in.)  in  inside  diameter  and  431.8  mm  (17° 
in-.  )"  in  length.  The  tube  was  heated  using  a  cartridge 
heater,  and  tLt^was  instrumented  with  8  thermocouples  to 
measure  the  wall  temperature.  A  Hewlett-Packard  3497A  data 
acquisition/ control  unit  and  a  9826A  computer  were  used  to 
collect  and  process  data.  The  condenser  was  cooled  by  an 
ethylene  glycol-water  mixture,  which  was  maintained  at  about 
-17  by  means  of  an  R- 12  refrigeration  system.  Nine  data 
runs  were  completed  to  de-bug  the  experimental  apparatus  and 
to  check  for  reproducibility.  During  all  data  runs,  espe¬ 
cially  at  higher  heat  fluxes  (greater  than  10  kW/ 
temperature  variations  were  observed  along  and 
active  boiling  length  of  the  test  tube.  The 
compared  with  data  found  in  the  literature  and'  reasonable 
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graduated  cylinder  is  supplied  with  oil  by  reservoir  4.  The 
valves  V5  ,  V6  and  V7  are  provided  for  emptying  the  evapo¬ 
rator  prior  to  dismantling  it,  for  example,  for  changing  the 
evaporator  tube.  Also,  fresh  R-114  can  be  added  to  the 
reservoir  while  valve  V5  is  closed  and  V6,  V7  are  open  and 
by  providing  R-114  vapor  through  V8  and  V10. 

The  test  tube  was  placed  horizontally  in  the  evaporator 
after  a  240  VAC  electrical  cartridge  heater  of  203.2  mm  (8 
in.)  in  length  had  been  inserted  in  its  middle  section.  The 
whole  system  was  connected  to  the  vacuum  pump  through  V8  and 
V10  in  order  to  remove  noncondensible  gases.  The  relief 
valve  was  set  to  a  gage  pressure  limit  of  138  kN/m2  (20 
psi),  which  is  50%  less  than  the  manufacturer- recommended 
working  pressure  limit  of  the  Pyrex  glass  tees.  All  connec¬ 
tions  were  assembled  with  Swagelok  fittings  sealed  with 
Teflon  ferrules  to  ensure  leak  tightness.  Also,  all  valves 
used  in  the  system,  except  control  valve  VC  and  by-pass 
valve  V9 ,  were  made  of  stainless  steel.  Copper  tubes  of 
9.5-mm  (3/8  in.)  diameter  were  used  for  all  piping,  except 
in  the  oil  section,  where  6.35-mm  (1/4  in.)  tubes  were  used. 


B.  BOILING  TEST  SECTION 


1.  Evaporator 


The  evaporator  is  a  T-shaped  container  made  of 
Corning  Pyrex  glass.  Figure  2.3  shows  a  schematic  of  the 
evaporator  with  dimensions.  The  selection  of  Pyrex  glass 
offers  several  advantages:  it  is  corrosion- resistant  and 
transparent  and  has  a  smooth  interior  surface  (this  mini¬ 
mized  nucleate  boiling  at  the  inner  surface  of  the 
container)  compared  to  metals,  and  is  more  stable  with 
temperature  and  pressure  variations  compared  to  ordinary 
glass.  The  operating  pressure  limit  of  the  evaporator  is 
guaranteed  by  the  manufacturer  up  to  a  gage  pressure  of  200 


II.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


A.  OVERVIEW  OF  THE  SYSTEM 

An  overall  schematic  representation  of  the  experimental 
apparatus  is  shown  in  Figure  2.1,  and  a  photograph  is  shown 
in  Figure  2.2.  The  apparatus  consists  of  a  Pyrex-glass 
evaporator  with  the  boiling  test  tube  1,  a  Pyrex-glass 
condenser  2,  an  R-114  liquid  reservoir  3,  an  oil  reservoir 
4,  a  graduated  oil  cylinder  5,  a  vacuum  pump  6,  a  cooling 
section  consisting  of  a  1/2-Ton  R-12  refrigeration  unit,  a 
tank  filled  with  water- ethylene  glycol  mixture  7,  and  a  pump 
8.  The  evaporator  1,  condenser  2,  oil-reservoir  4, 
graduated- oil  cylinder  5  and  R-114  reservoir  3,  are  placed 
in  a  Plexiglas  chamber  9 ,  where  the  temperature  can  be  main¬ 
tained  at  lower  values  (up  to  10  °C)  relative  to  ambient 
temperature.  This  chamber  also  served  as  a  safety  barrier  in 
case  of  an  inadvertent  overpressurization  of  the  glass 
apparatus . 

The  vapor  rising  from  the  evaporator  is  condensed  in  the 
condenser  2  and  the  liquid  is  then  fed  back  into  the  evapo¬ 
rator  through  valve  V5 .  In  this  manner,  the  specific 
enthalpy  of  condensation  ("latent  heat”)  of  R-114  vapor  is 
removed  in  the  condenser.  The  evaporator  of  the  R-12 
refrigeration  system  is  made  of  a  copper  coil,  which  is 
immersed  in  the  water- ethylene  glycol  tank  in  order  to  main¬ 
tain  the  temperature  of  the  mixture  at  about  -17  °C.  An 
8-GPM,  turbine-type  pump  8  transports  the  water-ethylene 
glycol  mixture  through  the  control  valve  VC  to  the 
condenser. 

The  oil  can  be  added  into  the  evaporator  1  using  the 
graduated  cylinder  5  by  opening  valves  VI  and  V4 .  The 


fluxes  and  mass  fractions.  They  observed  an  increase  in 
heat-transfer  coefficient  with  less  than  5%  mass  fraction  of 
oil  only  at  high  heat  fluxes  (such  as  higher  than  20  kW/m2 ) . 

The  effects  of  oil  on  the  heat  -  transfer  performance  of 
refrigerants  in  a  tube  bundle  have  been  summarized  by  Arai , 
et  al.  [Ref.  18]  in  the  following  manner:  "The  presence  of 
oil  in  a  shell-and- tube  evaporator  is  believed  to  have  two 
opposing  effects  on  the  heat  -  transfer  coefficient  of  the 
evaporator  tube  bundles.  Firstly,  many  people  suspect  that 
the  oil  left  out  from  the  evaporating  refrigerant  tends  to 
choke  the  surface  structure,  thereby  lowering  the  heat- 
transfer  coefficient.  Secondly,  from  the  experience  of 
conventional  heat  exchangers,  it  has  long  been  recognized 
that  the  presence  of  oil  enhances  foaming  on  the  upper  part 
of  the  tube  bundles.  This  foaming  increases  the  vigor  of 
refrigerant  movement  around  the  tubes,  so  that  heat  transfer 
is  enhanced." 

F.  THESIS  OBJECTIVES 

Based  on  the  foregoing  discussion,  it  is  clear  that 
little  information  exists  in  the  literature  to  enable  the 
successful  design  of  R-114  evaporators,  especially  when 
using  enhanced  boiling  surfaces,  in  the  presence  of  lubri¬ 
cating  oils.  The  objectives  of  this  research  effort  were 
therefore  to: 

1.  Design,  construct  and  instrument  an  experimental 
apparatus  in  order  to  investigate  the  heat  -  transfer 
performance  of  R-114  and  R-114-oil  mixtures  from 
enhanced  evaporator  tube  surfaces , 

2.  Successfully  operate  the  system  to  produce  repeatable 
data,  and 

3.  Obtain  smooth-tube  baseline  data  for  R-114  to  be  used 
as  a  standard  for  the  future  data  on  advanced 
surfaces . 
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boiling  heat-transfer  coefficients  decrease  with  increasing 
oil  concentration  in  R-ll,  R-12,  R-22,  R-113  and  R-114.  For 
oil  concentrations  less  than  3%,  however,  Dougherty  and 
Sauer  [Ref.  14]  demonstrated  that  the  heat-tansfer  coeffi¬ 
cient  may  be  increased  slightly  for  R-ll-oil  and  R-113-oil 
mixtures.  They  attributed  this  trend  to  an  increase  of 
foaming  of  the  mixtures.  This  result  also  has  been  confirmed 
by  Sauer,  Gibson  and  Chongrungreong  [Ref.  15]  using  R-12-oil 
mixtures,  and  by  Henrici  [Ref.  10]  using  R-114-oil  mixtures. 
However,  in  a  recent  experimental  study  using  R-113-oil 
mixtures,  Jensen  and  Jackman  [Ref.  13]  did  not  observe  such 
an  enhancement  with  low  oil  concentrations.  The  model  of 
the  governing  processes  involved  in  boiling  of  refrigerant- 
oil  mixtures  has  been  given  by  Stephan  [Ref.  16]  as  follows: 
"Obviously  the  reduction  in  heat  transfer  in  a  mixture  is 
due  to  the  fact  that  the  liquid  near  the  vapour- liquid 
interface  of  growing  bubbles  becomes  denuded  in  the  more 
volatile  components.  Thus,  a  concentration  gradient  exists 
between  the  interface  and  the  bulk  of  the  liquid,  and  the 
more  volatile  components  have  to  diffuse  through  the  liquid 
in  order  to  reach  the  interface.  A  mass- transfer  resistance 
thus  is  introduced,  which  is  not  present  when  boiling  single 
components . " 

The  effect  of  the  oil  viscosity  on  heat  transfer  in 
boiling  of  ref rigerant - oil  mixtures  has  been  considered  by 
Sauer,  Gibson,  and  Chongrungreong  [Ref.  15].  Oil  viscosity 
is  believed  to  significantly  contribute  to  the  changes  in 
boiling  characteristics  of  ref rigerant - oi 1  mixtures.  The 
heat-transfer  coefficient  is  greatest  for  the  highest 
viscosity  of  oil  for  the  same  temperature  difference. 

The  effect  of  oil  on  boiling  of  R-12  from  a  Gewa-T  tube 
was  investigated  by  Stephan  and  Mitrovic  [Ref.  17].  They 
reported  that  the  presence  of  oil  in  R-12  reduced  the  heat- 
transfer  coefficient  almost  over  the  entire  range  of  heat 
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refrigerant,  boiling  pressure,  tube  diameter,  surface  condi¬ 
tion  of  the  tube,  and  maybe  the  "submergence  effect" 
(resulting  from  the  hydrostatic  liquid  head)  on  the  evapora¬ 
tion  process.  Experimental  observations  have  shown  that 
changes  in  magnitude  of  these  properties  and  conditions 
significantly  affect  pool-boiling  heat-transfer ' performance . 
[Ref.  7]. 

E.  NUCLEATE  BOILING  PERFORMANCE  OF  REFRIGERANT -OIL  MIXTURES 

Most  of  the  refrigeration  systems  work  with  oil- 
lubricated,  hermetically  sealed  compressors,  which  are 
cooled  by  the  working  fluid.  Therefore,  the  fluids  being 
circulated  within  these  refrigeration  systems  generally  are 
not  pure  refrigerants,  but  are  refrigerant-oil  mixtures. 
The  saturation  temperature  of  a  refrigerant-oil  mixture  is 
higher  than  that  for  the  pure  refrigerant  at  the  same  pres¬ 
sure.  Also,  both  viscosity  and  surface  tension  increase 
when  oil  is  added  to  refrigerants.  Other  thermodynamic  prop¬ 
erties  of  refrigerants,  such  as  density  and  specific  heat, 
are  also  significantly  changed  by  the  addition  of  oil.  Thus, 
adequate  knowledge  on  the  influence  of  oil  on  the  boiling  of 
refrigerants  is  very  important  for  the  successful  design  of 
refrigeration  systems.  The  literature  reveals  only  very 
little  information  on  the  boiling  performance  of  R-114-oil 
mixtures  [Ref.  10].  In  particular,  heat-transfer  perform¬ 
ance  of  the  enhanced  evaporator  surfaces  in  refrigerant-oil 
mixtures  is  a  poorly-understood  subject  area.  Since  these 
advanced  boiling  surfaces  contain  a  high  density  of  pores, 
the  presence  of  oil  can  generally  result  in  a  greater 
performance  degradation  than  on  more- conventional ,  smooth 
surfaces. 

Various  studies  which  are  surveyed  and  listed  by  Jensen 
and  Jackman  [Ref.  13]  have  shown  that,  in  general,  nucleate 


transfer  coefficient  of  R-114  increases  remarkably  at  high 
pressure  near  the  critical  point.  Their  results  showed  that 
the  effect  of  surface  roughness  changes  with  pressure  and 
the  effect  gradually  vanishes  as  pressure  approaches  the 
critical  point.  Using  a  horizontal  smooth  copper  tube, 
Henrici  [Ref.  10]  investigated  boiling  heat-transfer 
performance  of  R-114  over  a  heat -flux  range  from  0.1  to  100 
kW/m2  and  over  a  pressure  range  from  37  to  252  kN/m2 . 

Stephan  and  Mitrovic  [Ref.  11]  investigated  the  perform¬ 
ance  of  pure  R-114  from  a  Gewa-T  enhanced  surface.  Boiling 
occurred  outside  a  horizontal  bundle  of  these  tubes.  They 
reported  that  the  position  of  the  tube  within  the  bundle 
does  not  effect  the  heat-transfer  coeff icient •  considerably . 
The  reason  for  this  effect  was  explained  in  the  following 
manner:  "The  flow  pattern  around  a  tube  seems  to  be  of  very 
weak  influence  and  the  heat-transfer  process  is  mainly 
governed  by  the  bubble  formation  and  the  two-phase  flow 
inside  the  T-shaped  channel.  The  top  part  of  the  fins,  which 
were  rolled  to  form  the  T,  practically  did  not  directly 
contribute  to  the  heat- transfer  process  because,  even  at 
maximum  heat  fluxes,  vapour  bubbles  were  not  produced 
there . " 

Marto  and  Hernandez  [Ref.  12]  also  investigated  nucleate 
pool-boiling  characteristics  of  a  Gewa-T  surface  in  R-113. 
Their  experimental  results  showed  a  disagreement  with 
Stephan  and  Mitrovic 's  above  explanation.  They  showed  that 
"the  liquid-vapor  motion  may  not  always  be  confined  to  flow 
circumferentially  within  the  channels  formed  by  the  T-shaped 
fins  as  previously  postulated.  Surface  imperfections  on  the 
tips  of  the  Gewa-T  fins  may  cause  bubble  nucleation  from 
these  sites.  The  presence  of  these  active  sites  may  further 
enhance  the  performance  of  the  Gewa-T  surface." 

The  factors  which  appear  to  affect  boiling  heat-transfer 
performance  of  a  single  tube  are  the  physical  properties  of 
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The  pool  boiling  heat-transfer  performance  of  R-113  from 
the  above-mentioned  commercially-available  enhanced  surfaces 
has  been  investigated  by  Marto  and  Lepere  [Ref.  6].  They 
reported  that  for  R-113,  the  enhanced  surfaces  showed  a  two 
to  tenfold  increase  in  the  heat-transfer  coefficient  when 
compared  to  a  plain  tube.  The  High-Flux  surface  was  most 
effective  over  a  broad  range  of  heat  fluxes,  whereas  the 
Gewa-T  surface  was  only  effective  at  high  heat  fluxes  (near 
100  kW/m* ) .  The  Thermoexcel-E  surface  showed  similar  gains 
in  the  heat-transf er  coefficient  to  that  of  the  High-Flux 
surface  below  10  kW/m* .  They  also  reported  that  the  degree 
of  superheat  required  to  activate  the  enhanced  surfaces  was 
less  than  the  plain  tube,  and  was  sensitive  to  initial 
surface  conditioning  and  fluid  properties. 

D.  NUCLEATE  POOL  BOILING  PERFORMANCE  OF  PURE  R-114 

The  heat- transfer  performance  of  pure  R-114  in  the 
evaporator  section  of  an  air-conditioning  unit  can  be  esti¬ 
mated  using  one  of  the  suggested  commonly-used  correlations 
for  nucleate  pool  boiling,  which  are  listed  by 
Chongrungreong  and  Sauer  [Ref.  7] .  However,  because  of  the 
lack  of  complete  understanding  of  the  nucleate  boiling 
mechanism,  these  correlations  should  be  compared  with  exper¬ 
imental  data  if  it  is  possible.  So  far,  the  literature 
reveals  a  limited  amount  of  nucleate  boiling  experimental 
investigations  using  R-114.  Happel  [Ref.  8],  using  a  hori¬ 
zontal,  pure-nickel  tube  with  a  surface  roughness  of  0.43 
micrometer,  investigated  the  boiling  characteristics  of  pure 
R-114  at  various  pressures.  The  effect  of  pressure  on 
nucleate  boiling  heat  transfer  in  R-114  was  also  investi¬ 
gated  experimentally  by  Nishikawa  and  his  co-workers 
[Ref.  9].  They  operated  over  a  reduced  pressure  (P/Pc  ) 
range  from  0.03  to  0.98,  and  demonstrated  that  the  heat 


a 


on  the  above  discussion,  R-114  is  probably  the  most  suitable 
refrigerant  for  use  on  board  submarines  and  ships. 

Further,  R-114  is  considerably  more  stable  toward  heat 
than  R-12.  All  standard  ferrous  and  nonferrous  metals, 
except  zinc  and  magnesium,  and  aluminum  alloys  containing 
appreciable  amounts  of  zinc  and  magnesium,  may  be  used  with 
R-114  [Ref.  3].  The  physical  properties  of  R-114  and  its 
comparison  with  the  other  common  refrigerants  are  summarized 
in  Appendix  A,  pressure-enthalpy  diagram  of  R-114  is 
provided  in  Appendix  B  and  the  range  of  application  of 
refrigerants  are  provided  in  Appendix  C  [Ref.  4]. 

C.  ENHANCED  BOILING  TUBE  SURFACES 

The  evalution  of  enhanced  surface  geometries  that 
promote  high-performance  nucleate  boiling  has  been  surveyed 
by  Webb  [Ref.  5].  Two  basic  types  of  enhanced  boiling 
surfaces  are  employed  on  the  outer  surface  of  tubes:  porous 
coatings  and  reentrant  grooves.  The  High-Flux  surface  is  a 
sintered,  porous  metallic  matrix  bonded  to  the  base  surface. 
The  Gewa-T  and  Thermoexcel-E  surfaces  have  reentrant  -  type 
nucleation  sites  formed  by  cold  working  an  integral  finned 
tube  wall.  See  Figures  1.1,  1.2  and  1.3  for  details. 

The  key  to  the  high  performance  of  these  surfaces  is 
attributed  to  three  factors  [Ref.  5]  :  "(1)  a  pore  or  reent¬ 
rant  cavity  within  a  critical  size  range,  (2)  interconnected 
cavities,  and  (3)  nucleation  sites  of  a  reentrant  shape. 
When  the  cavities  are  interconnected,  one  active  cavity  can 
activate  adjacent  cavities.  It  appears  that  the  dominant 
fraction  of  the  vaporization  occurs  at  a  very  thin  liquid 
film  within  the  subsurface  structure.  The  reentrant  cavity 
shape  provides  a  stable  vapor  trap,  which  will  remain  active 
at  very  low  liquid  superheat  values." 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  order  to  reduce  the  size  and  weight  of  air- 
conditioning  systems  and  water  chillers,  the  United  States 
Navy  recently  has  considered  two  different  approaches:  1)  to 
change  the  type  of  the  refrigerant  from  R-ll  to  R-114,  and 
2)  to  use  enhanced  evaporator  tube  surfaces,  such  as  Gewa-T, 
Thermoexcel-E ,  and  High-Flux.  Schematic  drawings  and  photo¬ 
micrographs  of  these  enhanced  surfaces  at  various  magnifica¬ 
tions  are  given  in  Figures  1.1  ,  1.2  and  1.3.  A  combination 
of  R-114  and  enhanced  surfaces  may  lead  to  significant 
improvements  in  system  efficiency  as  well  as  to  smaller  and 
lighter  components  in  future  Navy  applications. 

B.  ADVANTAGES  OF  USING  R-114 

The  latent  heat  per  unit  volume  of  R-114  is  50  BTU/ft* 
at  70  °F.  As  a  result,  the  heat  carried  by  each  cubic  foot 
of  R-114  is  nearly  twice  as  much  as  that  carried  by  R-ll 
[Ref.  1}  at  the  same  temperature.  ,  In  other  words,  the 
volume  of  R-114  vapor  circulated  per  ton  of  refrigeration  is 
about  half  that  for  R-ll.  Inherently,  the  size  of  the  system 
components  can  be  reduced  significantly  using  R-114  rather 
than  R-ll.  R-12  and  R-22  in  fact  can  transfer  the  heat  more 
efficiently  than  R-114,  but  they  require  higher  pressure 
ranges.  However,  operating  pressures  for  R-114  are  moderate, 
similar  to  R-ll,  such  that  light-weight  components  may  be 
desirable  and  less-experienced  personnel  may  be  involved  in 
the  installation  and  maintenance.  Whi*e  the  refrigerants  in 
general  are  nonflammable,  R-114  and  R-12  are  rated  to  have 
the  lowest  toxicity  among  all  refrigerants  [Ref.  2].  Based 
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temperature 

average  wal~  temperature  at  the  thermocouple 
location 

calibration  bath  temperature 
temperature  of  the  thermocouple  location 
saturation  temperature  of  the  boiling  liquid 
outer  wall  temperature  of  the  boiling  test  tube 
voltage  across  the  cartridge  heater 
voltage  reading  by  AC-DC  true  RMS  converter 
thermal  diffusivity 

volumetric  thermal  expansion  coefficient 
uncertainty  in  measurement  and  calibration 
viscosity 

kinematic  viscosity 

volume  fraction  of  pure  refrigerant 
saturated  liquid  density 
saturated  vapor  density 
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area 

tube  outside'  surface  area  of  active  boiling  section 
cross-sectional  area  of  tube 

tube  outside  surface  area  of  non-boilng  section 
conversion  coefficients  for  thermocouple 
specific  heat 
diameter 

tube  inside  diameter 
tube  outside  diameter 

diameter  at  the  position  of  the  thermocouple 

outer  diameter  of  the  boiling  tube 

discrepancy 

thermocouple  reading 

gravitational  acceleration 

heat- transfer  coefficient 

latent  heat  of  vaporization 

current 

current  reading  by  AC  Current  Sensor 
thermal  conductivity  of  liquid 
thermal  conductivity  of  copper 
active  boiling  tube  length 
non-boiling  length  of  the  test  tube 
Nusselt  number 
pressure 

critical  pressure 

tube  outside  wall  perimeter 

Prandtl  number 

heat- transfer  rate  from  boiling  surface 
heat-transfer  rate  through  one  unboiling  end 
heat -trans f er  rate  from  cartridge  heater 
Rayleigh  number 
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kN/m1  (30  psi).  Each  end  of  the  Pyrex-glass  evaporator  was 
fitted  with  a  cast-iron  flange  and  a  gasket.  A  detailed 
sketch  of  the  cast-iron  flange  is  shown  in  Figure  2.3.  Two 
aluminum  flanges,  210  mm  in  diameter  and  12.7  mm  in  thick¬ 
ness,  were  bolted  to  the  cast-iron  flanges.  Thus,  all 
fittings  were  connected  through  the  aluminum  flanges  to  the 
Pyrex-glass  evaporator. 

2.  Test  Section 


A  schematic  drawing  of  the  test  section  is  shown  in 
Figure  2.4.  The  boiling  test  tube  was  held  in  place  by  two 
Teflon  bushings,  which  were  attached  to  the  aluminum  flanges 
at  both  ends  of  the  T-shaped  evaporator.  Four  studs  were 
used  in  order  to  attach  the  Teflon  bushing  on  the  aluminum 
flange.  The  Teflon  bushing  and  both  ends  of  the  test  tube 
were  sealed  by  means  of  two  rubber  0-rings  in  order  to 
achieve  leak  tightness.  An  additional  0-ring  was  placed 
between  the  aluminum  flange  and  the  Teflon  plug  at  each  end. 

3 .  Boiling  Tube 

A  smooth,  hard-copper  tube,  15.9  mm  (5/8  in.)  in 
outer  diameter,  12.7  mm  (1/2  in.)  in  inside  diameter  and 
431.8  mm  (17  in.)  in  length,  was  used  to  provide  baseline 
data  for  the  nucleate  pool  boiling  heat-transfer  coefficient 
of  pure  R-114.  The  heater  was  a  1000-Watt  240-Volt 
stainless  -  steel  cartridge,  6.35  mm  (1/4  in.)  in  outer  diam¬ 
eter  and  203.2  mm  (8  in.)  in  length.  The  cartridge  heater 
was  inserted  into  a  copper  sleeve,  which  was  6.35  mm  (1/4 
in.)  in  inside  diameter,  12.7  mm  (1/2  in.)  in  outside  diam¬ 
eter  and  203.2  mm  (8  in.)  in  length.  In  order  to  provide  a 
uniform  heat  flux,  the  cartridge  heater  and  the  copper 
sleeve  were  soldered  together.  Pulido  [Ref.  19]  showed  that 
(using  a  similar  type  of  heater-element  construction)  the 
heater  element  provides  a  uniform  heat  flux  along  the 


circumferential  direction.  The  heater  element  used  during 
this  thesis  was  assumed  to  have  similar  characteristics . 
The  cartridge  heater  and  the  copper  sleeve  were  then 
inserted  as  a  unit  into  the  middle  portion  of  the  test  tube 
using  a  tight  mechanical  fit.  (Note:  A  shrink- fit  process 
is  essential  to  minimize  the  thermal  contact  resistance. 
However,  owing  to  the  time  constraint  in  this  thesis,  the 
shrink-fit  process  was  delayed  and  a  slide  fit  with  a  clear¬ 
ance  of  less  than  0.01  mm  was  used.  Thus,  the  resulting  tube 
was  used  mainly  for  the  purpose  of  obtaining  a  successfully- 
operating  apparatus).  The  active  boiling  length  of  the  test 
tube  was  203.2  mm  (8  in.)  in  the  middle  portion  of  the 
boiling  tube.  In  order  to  compute  the  actual  average  heat 
flux  in  the  heated  portion  of  the  tube,  a  suitable  correc¬ 
tion  was  applied  for  the  natural  convection  created  at  both 
ends.  See  Heat  Flux  Calculation  section  in  Chapter  4  for 
details . 

To  measure  the  wall  surface  temperature  of  the 
boiling  tube,  8  thermocouples  were  inserted  into  8  grooves 
which  were  machined  on  the  outside  of  the  copper  sleeve.  As 
shown  in  Figure  2.5,  these  thermocouples  were  located  at 
different  axial  and  circumferential  locations.  Four  thermo¬ 
couples  (1,2, 5, 6)  were  placed  with  90-degree  separation 
around  the  copper  sleeve  to  measure  the  temperature  distri¬ 
bution  at  the  mid  cross  section  of  the  active  length  of 
boiling  tube.  The  thermocouple's  wire  was  glued  down  at 
several  places  along  the  thermocouple's  channel  using  Epoxy. 
The  longitudinal  temperature  distribution  on  the  active 
boiling  section  was  also  measured  by  the  symmetrically- 
located  4  thermocouples  (3, 4, 7, 8).  The  exact  locations  of 
the  thermocouples  and  dimensions  of  the  thermocouple  grooves 
are  given  in  Figure  2.5.  All  thermocouple  grooves  were 
axially  machined  from  the  location  of  the  thermocouple  hot 
junctions  to  the  nearest  end  of  the  copper  sleeve. 


C.  CONDENSER  SECTION 


The  condenser  was  also  a  T-shaped  container  made  of 
Corning  Pyrex  glass.  It  was  identical  to  the  evaporator.  The 
position  of  the  condenser  can  be  seen  in  Figure  2.1.  R-114 
was  condensed  on  a  vertical  copper  coil,  which  was  inserted 
in  the  Pyrex-glass  condenser.  The  copper  coil  was  fabricated 
into  a  76.2  mm  (3  in.)  diameter  coil  using  9.5  mm  (3/8  in.) 
copper  tube.  The  active  condensation  length  was  estimated  to 
be  4.5  m  (15  f t ) . 

The  top  portion  of  the  condenser  was  connected  to  a 
portable,  mechanical  vacuum  pump  to  remove  noncondensible 
gases  from  the  apparatus.  The  bottom  of  the  condenser  was 
also  connected  to  the  evaporator  via  valve  V5  in  order  to 
return  the  condensed  R-114  liquid  to  the  evaporator.  The 
cooling  liquid,  i.e.,  water- ethylene  glycol  mixture,  entered 
the  top  portion  of  the  condenser,  through  the  copper  coil 
and  left  the  condenser  from  the  bottom  to  return  to  the 
water-ethylene  glycol  tank.  The  condenser  was  placed  verti¬ 
cally  and  connected  to  the  vapor  outlet  of  the  evaporator 
using  L-shaped  aliminum  tube,  50.8  mm  (2  in.)  in  diameter. 
The  maximum  vapor  velocity  of  R-114  vapor  through  the 
aluminum  tube  was  found  to  be  about  0.5  m/sec.  A  Bourdon 
gage  with  a  range  of  absolute  vacuum  to  a  gage  pressure  of 
1030  kN/m2  (150  psi)  and  a  relief  valve  which  was  set  to  138 
kN/m2  (20  psi)  were  placed  on  the  L-shaped  aluminum  tube. 

D.  OIL  ADDING  SECTION 

To  study  th».  boiling  performance  of  R- 114-oil  mixtures, 
a  cylindrical  aluminum  reservoir,  152.4  mm  (6. in)  in  diam¬ 
eter  and  152.4  cm  (6  in.)  in  height,  and  a  glass  oil 
cylinder  were  installed  above  the  evaporator.  The  relative 
positions  of  the  oil  reservoir  and  the  oil  cylinder  are 
shown  in  Figure  2.1.  The  oil  cylinder  was  355  mm  in  length 
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and  had  a  diameter  of  25.4  mm.  This  cylinder  was  specially 
ordered  to  achieve  a  resolution  of  0.5  ml.  The  oil  cylinder 
was  connected  to  the  oil  reservoir  through  valves  V3  and  V2 . 
The  addition  of  oil  into  the  evaporator  can  be  achieved 
through  VI  by  gravity  after  balancing  the  pressure  of  the 
oil  cylinder  with  that  in  the  evaporator  by  opening  valve 
V4 . 


E.  COOLING  SECTION 

1.  Water- Ethylene  Glycol  Mixture  Tank 

In  order  to  store  the  water- ethylene  glycol  mixture, 
a  special  tank  was  manufactured.  The  total  volume  of  the 
tank  was  0.154  m3  (0 . 48  m  x  0 . 48  m  x  0 . 66  m)  and  it  was  made 
of  12 . 7-mm- thick  Plexiglas  sheet.  All  sides  of  the  tank  were 
glued  together  with  methylene- chloride  solution.  The  joints 
were  held  together  with  small  screws  for  extra  strength.  The 
low  thermal  conductivity  of  Plexiglas  was  especially  suited 
to  minimize  heat  transfer  (from  room  to  water-ethylene 
glycol  mixture)  through  the  tank  walls.  The  tank  was  placed 
on  the  floor  and  all  sides  were  insulated  with  22  mm  (7/8 
in.)  thick  rubber  insulation  sheets.  The  cooling  mixture 
contained  49  liters  (13  gal.)  of  ethylene  glycol  and  94 
liters  of  (25  gal.)  distilled  water.  The  freezing  point  of 
this  mixture  was  about  -  25  °C. 

2 .  R- 12  Refrigeration  Plant 

A  1/2-Ton  R- 12  refrigeration  plant  was  installed  to 
cool  the  water- ethylene  glycol  mixture.  Figure  2  ..6  shows  a 
schematic  of  the  R- 12  refrigeration  plant.  It  consists  of  a 
compact-type  air-cooled  condenser,  a  compressor,  a  receiver, 
a  filter-drier  unit,  a  pressure  regulator,  a  pressure- 
control  switch  and  a  thermostatic  expansion  valve.  The 
evaporator  of  the  R-12  refrigeration  plant  was  constructed 


using  a  9.5  mm  (3/8  in.)  copper  tube,  which  was  immersed  in 
the  water-ethylene  glycol  tank.  The  temperature  of  the 
water-ethylene  glycol  mixture  is  controlled  by  both  a  ther¬ 
mostatic  expansion  valve  and  a  pressure  control  switch.  The 
R- 12  refrigeration  plant  was  adjusted  to  keep  the  tempera¬ 
ture  of  the  cooling  liquid  at  about  -17  °C. 

3.  Pump  and  Control  Valve 

An  8  GPM,  115  VAC  Burks  turbine- type ,  positive- 
displacement  pump  was  installed  on  the  floor  and  the 
suction  side  of  the  pump,  25.4  mm  (1  in.)  in  diameter,  was 
directly  coupled  to  the  water- ethylene  glycol  tank.  Cooling 
liquid  was  pumped  from  the  tank  to  the  condenser  through  the 
control  valve  VC.  Also,  a  by-pass  valve  V9  was  placed  before 
the  control  valve  VC  on  the  discharge  line.  The  use  of  the 
by-pass  line  served  two  important  purposes:  1,  it  avoided 
overloading  of  pump  8  in  the  event  value  VC  is  completely 
closed,  and  2,  it  provided  proper  mixing  for  the  "warm" 
stream  returning  from  the  condenser  2.  See  Figure  2.1  for 
the  positions  of  control  valve  VC  and  by-pass  valve  V9 .  The 
by-pass  valve  was  adjusted  and  set  permanently  so  that  a 
sufficient  range  of  flow  rates  can  be  achieved  through  valve 
VC  to  satisfy  proper  cooling  of  the  condenser. 

F.  R- 114  RESERVOIR 

An  aluminum  cylindrical  reservoir,  228  mm  (9  in. )  in 
diameter  and  254  mm  (10  in.)  in  height,  was  placed  verti¬ 
cally  between  the  evaporator  and  condenser  in  order  to  store 
R-114  as  a  liquid.  The  liquid  level  of  the  R-114  can  be 
observed  by  means  of  a  sight  hose  attached  on  the  reservoir 
with  proper  fittings.  R-114  reservoir  was  connected  to  the 
vapor  line  through  valve  V7  and  to  the  liquid  line  through 
valve  V6.  See  Figure  2.1  for  arrangement  of  the  reservoir. 
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An  aluminum  frame  (1.07  m  x  0.51m  x  0.61  m)  was 
constructed  to  locate  all  the  parts  of  the  apparatus,  except 
the  cooling  section.  All  four  vertical  sides  of  the  frame 
were  covered  with  12.7  mm  (1/2  in.)  thick  Plexiglas  sheets 
and  both  left  and  right  sides  were  provided  with  hinges  to 
enable  easy  access  to  the  components  of  the  apparatus. 
Aluminum  and  plywood  plates  were  used  to  cover  the  bottom 
and  top  sides  of  the  frame,  respectively.  The  valve  bodies 
of  VI  through  V8  were  placed  inside  of  the  front  Plexiglas 
sheet  with  the  valve  stems  penetrating  the  sheet.  Thus,  the 
valve  handles  were  accessible  from  outside  of  the  Plexiglas 
box.  The  whole  frame  was  placed  above  the  water-ethylene 
glycol  tank  with  aluminum  support  so  that  the  system  was 
very  compact . 

One  of  the  main  advantages  of  this  chamber  is  that  the 
temperature  surrounding  the  evaporator  can  be  reduced  rela¬ 
tive  to  the  ambient  temperature.  Also,  in  case  of  emergency, 
the  thick  Plexiglas  chamber  would  provide  a  safety  barrier 
to  personnel  and  equipment. 

H.  INSTRUMENTATION 

1 .  Power  Measurement 

A  240  Volt  AC  source  was  used  as  the  power  supply, 
and  it  was  adjusted  by  a  variac  in  the  range  of  0-260  Volt 
and  0-8  Ampere  according  to  the  desired  heat  flux  at  the 
surface  of  the  boiling  tube.  Power  input  to  the  boiling  tube 
was  measured  with  an  AC  current  sensor  and  an  AC-DC  true 
R.M.S  converter.  The  AC  current  sensor  was  connected  to  the 
input  line  of  the  heater  in  series  and  the  AC-DC  true  R.M.S 
converter  was  connected  in  parallel.  Figure  2.7  shows  a 
schematic  representation  of  the  power-measurement  devices. 
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Both  the  AC  current  sensor  and  the  AC-DC  true  R.M.S 
converter  were  connected  to  the  data  acquisition/reduction 
system . 


2 .  Temperature  Measurement 

Various  temperatures  were  monitored  throughout  the 
system  to  include: 

1.  Boiling  tube  wall  (8  thermocouples) 

2.  Liquid  temperature  (one  thermocouple) 

3.  Vapor  temperature  (two  thermocouples),  and 

4.  Water- ethylene  glycol  mixture  temperature  (one  ther¬ 
mocouple  ) 

The  locations  of  the  wall  thermocouples  in  the  sleeve  of  the 
boiling  tube  are  shown  in  Figure  2.5.  The  liquid  and  vapor 
thermocouples  were  inserted  into  the  two  specially- 
manufactured  thermocouple  wells.  Figure  2.8  shows  a  sche¬ 
matic  of  these  thermocouple  wells.  While  the  stainless  -  steel 
portion  minimizes  (owing  to  low  thermal  conductivity)  errors 
resulting  from  the  axial  conduction  of  heat  from  the 
surrounding,  the  copper  tip  helps  minimize  the  temperature 
drop  from  the  area  being  measured  to  the  thermocouple  loca¬ 
tion  (owing  to  the  high  thermal  conductivity  of  copper). 

All  the  temperature  measurements  were  accomplished 
by  0.245  mm  (30  gage)  copper-constantan  thermocouples.  Each 
thermocouple  measurement  was  read  directly  by  a 
Hewlett-Packard  3497A  data  acquisition  system,  which  was 
controlled  by  a  Hewlett-Packard  9826  computer.  Each  thermo¬ 
couple  was  scanned  for  0.8  seconds  and  twenty  readings  were 
averaged  to  obtain  a  more  accurate  measurement. 

A  total  of  five  thermocouples  were  calibrated.  Two 
thermocouples  were  made  from  the  beginning  of  a  spool  of 
copper-constantan  wire;  one  was  made  from  the  mid  portion 
and  two  were  from  the  end.  It  was  assumed  that  the  proper¬ 
ties  of  the  copper  and  constantan  do  not  change  along  a 
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a)  Corning  Pyrex  Glass  Evaporator  (D  x  Dj.  =  402x51  mm 
L  =  178  mm,  Lj  =  127  mm) 


b)  Cast  Iron  Flange  and  Gasket  (dj_=  190mm,  d0  =  210mm 
L  =  14mm ,  A  =  21°) 

Figure  2.3  Schematic  of  the  Pyrex  Glass  Evaporator. 


Schematic  of  the  Boiling  Test  Tube 


Figure  2.5  Positions  of  the  Thermocouples 
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and 


m2  =  h  p  /  kcAc 


(10) 


where : 


Tsat=  saturation  temperature  of  fluid  (°C) 

h  =  convective  heat - trans f er 
coefficient  (W/m2.K) 


Now,  assuming  negligible  heat  loss  from  the  thin  tube 
tip  (at  x  =  Lu) 


(d9  /  dx^{=I=i  0 


(ID 


and  also  at  the  base: 

Q  =  -k  A  (dT/dx)  =  -k,  A  (d0/dx)  (12) 

a  c  c  x  =  0  c  c  x  =  0 

whe  re : 

Qb  =  heat - trans fer  rate  from  the  base  (W) 

Using  the  above  boundary  conditions  (Eq.  11,  12),  the 

temperature  distribution  along  the  straight  fin  can  be 
expressed  as : 

(9/6b )  =  (  Cosh(m(Lu-  x))  /  Cosh(mLu)  )  (13) 
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location  (°C) 


Tav=  average  wall  temperature  at  the 
thermocouple  locations  (°C) 

1^  =  diameter  at  the  position  of  the 
thermocouple  (m) 

D2  =  outer  diameter  of  the  boiling  tube  (m) 

L  =  active  boiling  tube  length  (m) 

=  thermal  conductivity  of  the  copper  (W/m.K) 

1^0 =  outer  wall  temperature  of  the  boiling 
test  tube  (°C) 

Tb  =  temperature  at  the  base  of  straight  fin  (°C) 

The  temperature  distribution  of  a  fin  of  uniform  cross- 
sectional  area  is  given  by  [Ref.  20]  with  the  following 
assumptions: 

1.  One-dimensional  conduction  in  the  x  direction. 

2.  Steady-state  condition. 

3.  Constant  thermal  conductivity. 

4.  Negligible  radiation  from  the  surface  of  the  fin. 

3.  No  heat  generation  in  the  fin. 

6.  Uniform  convective  heat - 1 rans f er  coefficient  over  the 
fin  surface. 


d29  /  dx2  -  m2e  =  0 


(8) 


where : 


9(x)  =  T ( x )  -  Tsat 


(9) 
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where : 


(5) 


As  =  P  L 


A  =  cross-sectional  area  of  tube  (m2 ) 
c 

D  =  tube  outside  diameter  (m2 ) 
o 

D i  =  tube  inside  diameter  (m2 ) 

A  =  tube  outside  surface  area  of 
s 

nonboiling  section  (m2) 

p  =  tube  outside  wall  perimeter  (m) 

P  =  D  o 

Lu  =  non-boiling  length  of  the  test  tube  (m) 


It  was  assumed  that  the  temperature  at  the  base  of  the 
straight  fin  was  equal  to  the  average  wall  temperature  of 
the  active  boiling  section.  The  average  outer  wall  tempera¬ 
ture  of  the  active  boiling  section  was  calculated  using  the 
radial  conduction  equation  from  the  thermocouple  position  to 
the  surface  of  the  tube. 


vg 


8 

(  ^Tn)/8 


(6) 


Two=  Tavg  V  ln(P2/  ^)/(2  *  L  ]^)  )  (7) 


•T^  -  Two 


where : 


=  temperature  of  the  thermocouple 
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where : 


Vs  =  voltage  reading  by  AC-DC  true  RMS 
converter  (volts) 

Is  =  current  reading  by  AC  current 
sensor  (volts) 


The  total  length  of  the  boiling  test  tube  was  355.6  mm, 
while  the  active  boiling  length;  i.e.,  the'  length  of  the 
cartridge  heater,  was  203.2  mm.  The  geometry  of  the  boiling 
test  tube  is  shown  schematically  in  Figure  4.1. 
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Figure  4.1  Geometry  of  the  Boiling  Test  Tube. 

Natural- convection  heat  transfer  was  assumed  to  occur  at 
both  ends  of  the  test  tube  (the  distance  ) .  This  part  of 
the  tube  was  considered  to  be  a  straight  fin  of  uniform 
cross  section.  For  this  cylindrical  geometry: 


A  =  ir  (D*  -  D2)/4 
c  1 


(4) 
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11.  For  each  data  series,  seven  different  heat  fluxes  (1, 
2,  5,  10,  20,  50  and  95  kW/m2 )  were  selected. 

C.  HEAT- FLUX  CALCULATION 

According  to  the  description  of  the  boiling  tube 
construction,  the  cartridge  heater,  which  was  inserted  into 
the  boiling  tube  is  the  heat  source.  The  variations  in  heat 
flux  are  made  through  different  voltage  settings  of  the 
variac . 

%  =  VI  (1) 

where : 

Qh  =  heat- transfer  rate  from  the  cartridge 
heater  (W) 

V  =  Voltage  across  the  cartridge  heater 
element  (volts) 

I  =  current  through  the  heater  element  (amps) 

This  circuit  is  connected  in  series  with  an  AC  current 
sensor  and  in  parallel  with  an  AC-DC  true  RMS  converter. 
Each  of  these  two  sensors  produces  a  DC  output  in  the  range 
of  0-10  V,  which  is  automatically  read  by  the  data 
acquisition/ control  unit.  The  calibration  equations  of  these 
sensors  give  the  actual  values  as  follows: 


4.  Any  accumulated  noncondensible  gases  were  evacuated 
by  the  portable  mechanical  vacuum  pump  through  valve 
V8  and  Vll  (see  Figure  2.1). 

5.  The  iteractive  computer  program  DRPR1  was  operated 
and  the  desired  heat  flux  and  saturation  temperature 
were  given  to  the  program  as  reference  values. 

6.  The  cartridge-heater  voltage  was  then  adjusted  by 
variac  in  order  to  maintain  the  desired  heat  flux. 
The  desired  and  actual  heat  fluxes  were  compared 
continiously  by  the  computer  program  until  they 
agreed  to  within  2  percent. 

7.  The  actual  versus  desired  saturation  temperature  of 
the  liquid  R-114  was  also  monitored  by  computer 
program.  The  amount  of  cooling  liquid,  which  was 
being  circulated  through  the  copper  coil  in  the 
condenser,  was  regulated  by  control  valve  VC  in  order 
to  obtain  nearly  constant  saturation  temperature 
(±0.2  °C)  at  a  given  heat  flux. 

8.  For  all  the  consecutive  settings  during  both 
increasing  and  decreasing  heat  fluxes,  the  boiling 
was  allowed  to  stabilize  for  five  minutes  at  each 
power  setting  and  raw  data  (thermocouple  readings,  AC 
current  sensor  readings  and  AC-DC  true  RMS  converter 
readings)  were  recorded  in  a  user- specif ied  file. 

9.  At  a  given  power  setting  and  a  saturation  tempera¬ 

ture,  the  following  processed  data  were  recorded  as  a 
printout:  wall  temperatures  of  the  boiling  section, 

liquid  bulk  temperature,  vapor  temperature,  net  heat 
flux,  temperature  of  the  water- ethylene  glycol 
mixture,  wall  superheat  and  heat-transf er  coefficient 
of  the  R-114. 

10.  For  each  data  set,  the  above  procedure  beginning  with 
step  5  was  repeated. 
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position  (20  mm  above  the  boiling  test  tube)  with  the 
following  procedure  (See  also  Figure  2.1). 

1.  The  temperature  of  the  water- ethylene  glycol  tank  was 
reduced  to  about  -17  °C  by  the  R- 12  refrigeration 
unit . 

2.  The  apparatus  was  evacuated  to  about  29  in  Hg  vacuum 
by  the  vacuum  pump  (the  valves  Vll  and  V8  were  open 
and  V10  was  closed). 

3.  Upon  closing  Vll,  the  vapor  outlet  of  the  R-114 
supply  cylinder  was  connected  to  V10. 

4.  The  cooling  section  pump  8,  delivered  the  cooling 
liquid  from  the  water- ethylene  glycol  tank  to  the 
condenser  through  the  control  valve  VC. 

5.  The  R-114  vapor  condensed  on  the  condenser  coils  and 
the  liquid  was  collected  in  the  evaporator  by 
gravity,  while  V5  was  open  and  V6  was  closed. 

6.  Pump  8  was  then  stopped  when  the  desired  liquid  level 
of  the  evaporator  was  achieved. 

7.  The  pressure  of  the  apparatus  was  then  allowed  to 
increase  up  to  the  saturation  pressure  corresponding 
to  the  ambient  temperature. 

B.  NORMAL  OPERATION 

The  following  procedure  was  established  to  obtain  the 
heat- transfer  coefficient  of  R-114  from  the  smooth  copper 
test  tube: 

1.  The  test  tube  was  immersed  in  the  pool  (20  mm  below 
the  liquid  level). 

2.  The  R- 12  refrigeration  unit  was  operated  24  hours  in 
advance  in  order  to  reduce  the  temperature  of  the 
water- ethylene  glycol  tank  to  about  -17  °C. 

3.  The  data  acquisition/ control  unit,  computer,  power 
supply  and  cooling  section  pump  were  switched  on. 
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IV.  EXPERIMENTAL  PROCEDURE 


A.  PREPARATION 

1.  Pressure  Test  of  the  Apparatus  with  Air 

Upon  assembling  all  the  components,  including  the 
boiling  test  tube,  the  experimental  apparatus  was  pressur¬ 
ized  with  air  up  to  100  kN/m2  (15  psi).  A  soap-bubble  test 
was  first  carried  out  to  locate  any  leaks  from  each  compo¬ 
nent  of  the  apparatus.  All  detected  leaks  were  succesfully 
and  systematically  fixed  at  the  end  of  each  pressure  test. 

2.  Pressure  Test  of  the  Apparatus  with  R- 114 

A  second  pressure  test  was  performed  using  R-114 
vapor.  The  apparatus  were  charged  with  R-114  vapor  in  the 
following  manner.  The  apparatus  was  evacuated  to  25  in.  Hg. 
by  a  portable  mechanical  vacuum  pump  through  valves  Vll  and 
V8  (See  Figure  2.1  for  the  configuration  of  valves).  Upon 
closing  valve  Vll,  a  vapor  outlet  of  the  R-114  supply 
cylinder  was  connected  to  the  apparatus  by  means  of  valve 
V10.  Thus,  the  pressure  of  the  apparatus  was  increased  up  to 
the  saturation  pressure  of  R-114  at  that  ambient  temperature 
(ambient  temperature  was  21  °C  and  the  corresponding  satura¬ 
tion  pressure  (gage)  was  86  kN/m2  (12.5  psi)).  An  Automatic 
Halogen  Leak  Detector,  TIF  5000,  was  used  to  detect  R-114 
leakage  from  the  apparatus.  The  sensitivity  of  this  detector 
is  1/2  oz.  per  year,  3  PPM  max.  A  few  small  leaks  were 
observed  and  successfully  isolated. 

3.  Charging  the  Apparatus  with  R-114 

Upon  completing  the  pressure  tests,  the  evaporator 
of  the  apparatus  was  filled  with  R-114  liquid  up  to  a  marked 


Photograph  of  Data  Acquisition/Reduction  System. 


Appendix  E  shows  a  listing  of  computer  program  DRPR1  and 
subprogram  PLOT  and  Appendix  F  shows  an  example  of  represen¬ 
tative  data  run. 


TABLE  I 

HP  3497A  Channel  Assignments 


Channel 

Assignment 

25  - 

32 

Boiling  tube  wall  temperature 

33 

Liquid  temperature 

34  - 

35 

Vapor  temperature 

36 

Temperature  of  cooling  liquid 

62 

AC-DC  true  RSM  converter 

63 


AC  current  sensor 


STEPWISE  DATA- COLLECTION  AND  SOLUTION  PROCEDURE 


1.  Select  tube  type  (all  dimensions  of  the  boiling  test 
tube  are  included). 

2.  Set  desired  heat  flux  and  saturation  temperature  of 
the  boiling  liquid.  Wait  for  steady-state  conditions 
(See  Chapter  4  section  B  for  details). 

3.  Scan  all  channels  listed  in  Table  I  (thermocouples, 
AC  current  sensor  and  AC-DC  true  RMS  converter  read¬ 
ings  ) . 

4.  Save  these  raw  data  in  a  user- specif ied  file. 

5.  Convert  these  raw  data  readings  to  corresponding 
units  (temperature  in  °C,  current  in  Ampere,  voltage 
in  volt ) . 

6.  Compute  the  heat  -  transfer  rate  from  the  cartridge 
heater  (See  Chapter  4,  Section  C  for  details). 

7.  Compute  the  average  wall  temperature  of  the  boiling 
test  tube  and  calculate  the  wall  superheat  of  this 
data  set  (See  Chapter  4,  Section  C  for  details). 

8.  Compute  the  physical  properties  of  R-114  using  given 
correlations  at  film  temperature  (See  Appendix  H  for 
details ) . 

9.  Compute  the  natural- convection  heat  -  transfer  coeffi¬ 
cient  of  R-114  from  non-boiling  ends  of  the  test 
tube . 

10.  Compute  heat  losses  from  non-boiling  ends. 

11.  Calculate  the  heat  flux  from  boiling  test  tube  to  the 
boiling  liquid. 

12.  Calculate  boiling  heat-transfer  coefficient  of  the 
R-114  from  test  tube. 

13.  Store  the  heat  flux  versus  wall  superheat  values  for 
each  data  set  in  a  user- specif ied  plot  file. 

14.  Use  the  subprogram  PLOT  and  plot  the  data  run. 
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III.  DATA  ACQUI S ITION / REDUCTION 

A.  DATA  ACQUISITION  AND  STORAGE 

A  Hewlett-Packard  3497A  automatic  data  acquisition/ 
control  unit  was  used  to  read  temperatures  from  the  thermo¬ 
couples  and  to  read  current  and  voltage  values  of  the 
cartridge  heater  from  the  AC  current  sensor  and  the  AC-DC 
true  RMS  converter,  respectively.  A  Hewlett-Packard  9826A 
computer  was  used  to  control  the  data  acquisition/ control 
unit  and  to  analyze  and  store  data.  Figure  3.1  shows  a 
photograph  of  the  data  acquisition/ control  unit  and 
computer. 

Information  was  entered  through  the  keyboard  to  prompt 
the  data  acquisition/control  unit  to  automatically  scan  each 
channel.  Channel  assignments  are  listed  in  Table  I.  These 
raw  data  were  immediately  processed  and  a  hard-copy  printout 
was  provided.  Also,  these  data  were  transferred  to  a 
computer  disk  under  a  user- specif ied  file  name  to  keep  a 
permanent  record.  The  ability  to  store  raw  data  directly 
enabled  these  data  to  be  reduced  at  any  time  and  allowed 
flexibility  for  changes  to  data-reduction  software. 

B .  DATA  REDUCTION 

Following  data  acquisition  for  each  data  point,  results 
were  computed  according  to  the  stepwise  procedure  outlined 
in  the  next  section,  and  then  printed  on  a  Hewlett-Packard 
2671G  thermal  printer.  Heat  flu;:  versus  wall  superheat  (the 
temperature  of  the  boiling  surface  minus  the  liquid  tempera¬ 
ture)  were  also  stored  in  a  user- specif ied  plot  data  file 
for  subsequent  plotting  using  the  subprogram  PLOT.  A 
Hewlett-Packard  7470A  plotter  was  also  interfaced  with  the 
computer . 
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Sketch  of  a  Thermocouple  Well 


Figure  2.7  Schematic  of  the  Power  Measurement. 
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and  the  fin  heat-transf er  rate  is: 


h  p  kcAc  9bTanh(mLu) 


(14) 


where : 


Qf  =  heat - trans f er  rate  through  one  non-boiling 
length  of  the  test  tube  (W) 

The  average  difference  between  the  wall  temperature  and 
the  liquid  temperature  may  now  be  determined  from  the 
following  equation: 

fL“ 

(0/0b)  =  ( j(Cosh(m(Lu-  x))/Cosh(mLu)  )  dx  )  /  L u  (15) 

-'o 

0  =  (0b/  mLu )  Tanh(mLu5 


‘wo  A  sat 


where: 


T^o  "  ^at=  average  difference  between  outer 

wall  temperature  and  liquid  saturation 
temperature  (K) 

A  free  convection  correlation  stated  by  Churchill  and 
Chu  [Ref.  20]  was  applied  for  the  average  Nusselt  number  on 
a  horizontal  cylinder: 


0.387  Rar 


0. 60+- 


[1  +  (0.559/Pr)9/16]8/2?  I 


(16) 


The  average  Nusselt  number  is: 


where : 


h  D 


(17) 


k  =  thermal  conductivity  of  R-114  (W/m.K) 


Solving  for  h: 


h  - 


k  i 

0.60  + 


0.387  Ra*/6 
o 


. . I 

°o  I  [1  +  <0.559/Pr)9/16] 8/27  ) 


(18) 


Pr  =  — 

a 


(19) 


Ra, 


9BlfWn~TSAT)Do 


v  a 


(20) 


g  =  -  I  *  a 

p  LT 


(21) 


where : 


Pr  =  Prandtl  number 
v  =  kinematic  viscosity  (in2/s) 
a  =  thermal  diffusivity  (m2/s) 
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average  Rayleigh  number 


gravitational  acceleration  (m/s1) 


volumetric  thermal  expansion 
coefficient  (1/K). 


density  of  R-114  liquid  (kg/mJ) 


Now,  substitution  of  Equations  10,  15,  20  into  Equation 

18  results  in: 


gSD^P  Tanhir^-)  Lull  1/6 

’  o  B  |  >v'V  t 


0.60  +  0 . 387 


r  t  hP  ^ 1/2 
U  ^ 


(22) 


(1  +  (0. 559/Pr)9/16) 8/27 


Equation  22  is  solved  for  h  by  an  iterative  technique 
within  a  range  of  precision  of  0.001.  Knowing  the  value  of 
the  natural  convection  heat- transfer  coefficient  along  the 
non-boiling  ends,  the  total  heat-loss  rate  is  calculated 
from  Equation  14,  and  the  heat-transfer  rate  through  the 
boiling  section  is  obtained  by  subtracting  the  total  heat- 
loss  rate  from  the  total  heat- transfer  rate: 


^03%  2xQT 


(23) 


Q  =  V  2^f 


(24) 


V.  RESULTS  AND  DISCUSSION 
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A.  OUTLINE  OF  THE  DATA  RUNS 

Using  the  procedure  outlined  in  Chapter  3,  nine  data 
runs  were  completed,  primarily  to  de-bug  the  experimental 
apparatus  for  its  successful  operation.  For  all  data  runs, 
the  same  smooth  copper  tube,  described  in  Chapter  2,  was 
used  as  the  boiling  surface.  Each  data  run  consisted  of 
seven  different  heat  fluxes  (1000,  2000,  5000,  10000,  20000, 
50000  and  95000  W/m2)  with  a  specified  saturation  tempera¬ 
ture  (0,  10  or  20  °C) .  All  the  raw  data  were  stored  in  a 
user-specified  file,  named  "WHxx" ,  where  "W"  indicates 
Wieland  tube,  "H"  indicates  hard  copper  and  "xx"  indicates 
the  data  run  number.  Also,  the  heat-transfer  coefficient 
versus  heat  flux  data  were  stored  in  a  plot  file  named 
"Pxx" .  A  summary  of  these  data  runs  is  given  in  Table  II. 

B.  LONGITUDINAL  AND  CIRCUMFERENTIAL  TEMPERATURE  VARIATIONS 

During  all  data  runs,  especially  at  high  heat  fluxes 
(greater  than  10  kW/m2 ) ,  large  temperature  variations  were 
observed  along  the  active  boiling  length  of  the  test  tube. 
Figure  5.1  shows  the  measured  wall  temperatures  at  each  of 
the  thermocouple  locations,  while  Figure  5.2  shows  the 
temperature  variations  in  the  circumferential  direction  (at 
the  axial  mid  point)  at  two  different  heat  fluxes  (2  and  92 
kW/m2 ) .  Despite  the  uniform  heat  flux  provided  by  the 
heater,  the  temperature  measurements  showed  considerable 
variations  in  both  axial  and  circumferential  directions. 
The  reason  for  these  temperature  variations  may  be  explained 
as  follows:  As  explained  in  Chapter  2,  during  the  manufac¬ 
turing  of  the  boiling  test  tube,  the  copper  sleeve  was 
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inserted  into  the  boiling  test  tube  without  any  shrink-fit 
process.  The  slip-fit  process  (with  a  diameter  clearance  of 
about  0.01  mm)  used  in  this  study  was  believed  to  have 
caused  considerable  thermal  contact  resistance  at  the  inter¬ 
face  between  the  boiling  tube  and  the  copper  sleeve.  Also, 
it  is  quite  possible  that  this  contact  resistance  was  not 
uniform  along  both  the  axial  and  circumferential  directions. 
The  circumf erent ial  temperature  distribution  may  also  be 
attributed  to  a  circumferential  variation  in  boiling  heat- 
transfer  coefficient.  During  the  high  heat  flux  boiling 
regime,  it  was  observed  that  the  rate  of  bubble  formation  at 
the  top  of  the  horizontal  boiling  tube  was  lower  than  at  the 
rest  of  the  boiling  tube;  i.e.,  the  boiling  heat- transfer 
coefficient  at  top  of  the  boiling  tube  was  relatively  lower 
than  the  other  locations.  This  poorer  performance  at  the  top 
of  the  tube  may  be  attributed  to  the  fact  that  this  portion 
was  not  receiving  sufficient  amount  of  fresh  liquid  to 
replace  the  liquid  that  evaporated.  The  obstruction  of  the 
liquid  flow  could  be  caused  by  the  growing  vapor  blanket 
from  the  bottom  to  the  top  of  the  -tube.  Based  on  this  obser¬ 
vation,  it  is  expected  that  the  temperature  at  top  of  the 
boiling  tube  could  be  greater  than  the  other  locations.  As 
shown  in  Figure  5.2,  location  5,  which  represents  the  top  of 
the  horizontal  boiling  tube,  measured  the  highest  tempera¬ 
ture  difference  as  expected  from  the  above-mentioned  obser¬ 
vation. 

However,  positions  1  and  2  showed  a  sizable  discrepancy 
at  high  heat  flux,  and  this  was  due  most  likely  to  the  non- 
uniform  contact  resistance.  Also,  a  comparison  of  thermocou¬ 
ples  2,  7  and  8  show  a  similar  discrepancy.  It  is  clear 
therefore,  that  the  shrink- fit  process  must  be  used  in  the 
future . 


C.  PLOT  ANALYSIS  OF  NUCLEATE  BOILING  REGIME 


Figure  5.3  shows  a  typical  nucleate  pool-boiling 
performance  curve  of  the  smooth  copper  tube  in  R-114.  The 
observed  behavior  of  this  process  is  analyzed  by  studying 
the  different  regions  of  the  boiling  curve. 

From  point  A  to  point  B,  a  continuous  increase  in  wall 
superheat  (T^0  -lLat)  is  observed  when  heat  flux  is  increased. 
Also,  during  the  experimental  runs,  no  bubbles  were  observed 
along  the  boiling  tube  in  this  region  of  the  curve.  This 
region  corresponds  to  the  natural - convection  process. 

From  point  B  to  point  C,  a  reduction  in  wall  superheat 
is  observed,  while  the  heat  flux  continuously  increased. 
While  the  middle  portion  of  boiling  tube  showed  bubbles,  the 
rest  of  the  boiling  tube  showed  no  bubbles  during  the  exper¬ 
iment  in  this  region.  This  region  is  known  as  the  mixed 
boiling  region,  where  transition  from  natural  convection  to 
nucleate  pool  boiling  heat  transfer  takes  place.  The  wall 
superheat  continued  to  decrease  until  numerous  nucleation 
sites  became  active  (point  C). 

After  point  C,  the  wall  superheat  starts  to  increase 
with  increasing  heat  flux  as  shown  by  region  C  to  D.  A  very 
high  density  of  bubble  formation  was  observed  along  the 
active  boiling  section  of  the  test  tube  (only  a  few  bubbles 
were  observed  on  the  nonboiling  sections  of  the  test  tube). 
See  Figure  5.4.  This  observation  confirmed  the  validity  of 
the  natural  convection  model  assumed  for  data  reduction. 
However,  the  R-114  liquid  returning  to  the  evaporator  from 
the  condenser  caused  a  mild  turbulance  in  the  entrance 
region  of  this  return  line  in  the  evaporator.  Inherently, 
this  affects  the  natural- convective  heat- transfer  perform¬ 
ance  of  nonboiling  ends  of  the  test  tube. 

When  heat  flux  is  gradually  decreased,  the  curve  follows 
a  different  path  (from  D  to  E)  as  shown  in  Figure  5.3.  This 
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is  due  to  the  existing,  stable  nucleation  sites  remaining 
active  for  a  wide  range  of  heat  fluxes.  Similar  analysis  can 
also  explain  Figure  5.5,  which  represents  the  same  data  run 
but  on  a  different  basis  (heat  -  transfer  coefficient  versus 
heat  flux) . 

D.  REPRODUCIBILITY  TEST  OF  THE  APPARATUS 

In  order  to  test  the  reproducibility  of  the  experimental 
apparatus,  two  runs  were  performed  on  two  different  days 
(data  run  number  5  and  7  in  Table  II)  at  the  same  saturation 
temperature.  Figure  5.6  and  Figure  5.7  show  the  comparison 
of  these  two  data  runs .  It  can  be  seen  from  these  figures , 
that  there  is  very  good  agreement  between  these  two 
different  runs.  This  agreement  shows  the  ability  of  the 
apparatus  to  reproduce  data  runs  revealing  successful 
operation. 


E.  BOILING  PERFORMANCE  OF  SMOOTH  COPPER  TUBE  IN  R-114 

In  order  to  test  the  validity  of  the  data  taken  from  the 
present  experimental  apparatus  (despite  the  presence  of 
errors  owing  to  the  thermal  contact  resistance  mentioned 
earlier),  an  attempt  was  made  to  compare  the  present  data 
with  data  found  in  the  literature.  For  this  purpose,  data 
run  number  5  (10  °C  saturation  temperature  with  decreasing 

heat-flux  rate)  was  compared  with  two  sources  in  the  litera¬ 
ture  as  discussed  below. 

1.  Comparison  with  Chonerunereone- Sauer  Correlation 


d  by  Ch 
iling  p 
xtures 
s  based 
urces : 

231  . 


where : 


h  =  heat-transfer  coefficient  (W/m2.K) 

Q  =  heat  rate  (W) 

A  =  surface  area  (m2) 

D  =  characteris tic  dimension  of  heated  surface  (m2 ) 
UL  =  viscosity  of  saturated  liquid  (gr/m.s) 
hfg  =  latent  heat  of  vaporization  (W.s/g) 

Cp=  specific  heat  of  liquid  (J/g.K) 

kL=  thermal  conductivity  of  liquid  (W/m.K) 

P  =  pressure  (atmospheres) 

^  =  volume  fraction  of  pure  refrigerant 
^  =  liquid  density  (g/cmJ) 

Pv  =  vapor  density  (g/cmJ) 


Eq'n  (26)  agrees  very  well  with  theirdata  [Ref.  7] 
for  R-ll  as  well  as  with  Stephan's  data  [Ref.  22]  for  oil- 
free  refrigerants,  R-ll,  R-12,  R-13,  R-21,  R-22,  R-113  and 
R-114  with  errors  smaller  than  16  percent.  Figure  5.8  pres¬ 
ents  a  comparison  of  the  current  experimental  data  with  eq'n 
(26).  As  can  be  seen  from  this  figure,  the  current  data  are 
in  excellent  agreement  with  this  predictive  equation.  (The 
analysis  given  in  Appendix  I  shows  that  the  uncertainty  in 
the  boiling  heat - trans f er  coefficient  for  this  experimental 
investigation  is  about  2Q70.  )  This  agreement  was  somewhat 
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unexpected  since  the  present  data  includes  some  kind  of 
contact  resistance  in  it.  If  this  contact  resistance  were 
known,  and  were  removed  from  the  data,  a  higher  boiling 
coefficient  would  result. 

2.  Comparison  with  Data  of  Henrici 

The  nucleate  boiling  heat-transfer  coefficient  of 
R-114  and  R-114-oil  mixtures  from  a  horizontal,  smooth 
copper  tube  (30  mm  in  outer  diameter)  has  been  experimen¬ 
tally  investigated  by  Henrici  [Ref.  10].  Since  Henrici 's 
data  were  taken  from  a  tube  with  a  larger  diameter  than  the 
current  boiling  tube  diameter,  his  data  were  revised  to 
include  this  difference.  For  this  purpose,  a  semi- emprical 
correlation  to  show  the  effect  of  tube  diameter,  as  devel¬ 
oped  by  Cornwell  et  al.  [Ref.  21],  was  used: 

Nu  =  C  Re2 '3  (27) 


•where : 

Re  =  (q  D)  /  (hfgU  )  and  Nu  =  (h  D)  /  k 
C  =  150  for  refrigerant 


Using  this  relationship,  Henrici’ s  data  were  revised 
for  a  diameter  of  15.875  mm  (5/8  in.)  used  in  this 
experiment . 


h(D=15 .875  mm)  =  h(D=30  mm)  ( 30/ 15 . 875  ) 1  ' 3 
h(D= 15.875  mm)  =  1.24  h(D=30  mm) 
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A  line  representing  the  revised  Henrici's  data  is 
plotted  in  Figure  5.9  together  with  the  unrevised  data  line. 
It  can  be  seen  that  the  boiling  heat-transfer  coefficient 
increases  by  about  25  percent  when  the  diameter  is  decreased 
from  30  mm  to  15.875  mm  (50-percent  reduction  in  diameter). 

Figure  5.10  provides  Henrici's  natural  convection 
data  and  Henrici's  revised  boiling  data  in  comparison  with 
the  present  experimental  data.  As  expected,  the  natural- 
convective  heat - trans f er  coefficient  is  lower  than  the 
boiling  heat- transfer  coefficient.  The  boiling  heat-transfer 
coefficient  from  Henrici’s  data  is  higher  than  the  coeffi¬ 
cient  found  in  this  study.  The  reason  for  this  may  be 
explained  as  follows.  In  Henrici's  experimental  investiga¬ 
tion,  the  measurement  of  the  surface  temperature  has  been 
done  using  thermocouples  outside  of  the  boiling  surface. 
Unfortunately,  he  did  not  describe  in  detail  how  he  measured 
the  wall  temperatures.  If  the  thermocouples  were  not 
totally  imbedded  in  the  wall,  they  may  produce  incorrect 
wall  temperature  values.  If  a  portion  of  a  thermocouple  is 
in  direct  contact  with  the  liquid,  this  thermocouple  will 
measure  a  value  lower  than  the  actual  wall  temperature  and 
such  an  error  would  decrease  the  wall  superheat,  thus 
resulting  in  a  higher  heat- transfer  coefficient.  Also,  if 
the  contact  resistance  effect  is  removed  from  the  present 
data,  then  the  heat  transfer  coefficient  would  increase 
toward  the  Henrici  data. 

F.  EFFECT  OF  PRESSURE 

The  effect  of  pressure  on  the  nucleate-boiling  heat- 
transfer  performance  of  R-114  was  investigated  for  reduced 
pressures  (P/Pc  )  of  0.0268,  0.0391  and  0.0556  (data  run 
numbers  5,  8  and  9,  respectively  in  Table  II).  Figure  5.11 
shows  the  relative  increase  in  boiling  heat  -  transfer 
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coefficient  with  increase  of  pressure  as  expected.  In  order 
to  calculate  the  pressure  effect  on  the  boiling  performance, 
a  simple  equation  developed  by  [Ref.  7]  was  used. 

h  =  6.17  (q)0-55  ( (j))3-ss  P0'24  (28) 

where : 

h  =  heat- transfer  coefficient  (W/m2.K) 
q  =  heat  flux  (W/m2 ) 

0  =  volume  fraction  of  pure  refrigerant 
P  =  boiling  pressure  (atmospheres) 

Using  this  equation,  the  pressure  effect  on  the  boiling 
heat  transfer  coefficient  can  be  written  as  follows: 


h(P[)  =  h(P2)  (Px  /P2  )0-24 


According  to  the  above  equation,  the  effect  of  pressure 
from  127  kN/m2  to  87  kN/m2  decreases  the  boiling  heat- 
transfer  coefficient  by  about  ten  percent,  which  agrees 
reasonably  well  with  the  current  experimental  data. 
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TABLE  II 

Summary  of  the  Data  Runs 
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APPENDIX  A 


PHYSICAL  PROPERTIES  OF  FREON  FLUOROCARBON  COMPOUNDS 


"FREON**  11 

"FREOfT  U 

"FREON'*  113 

"FREON"  114 

Cbcmieal  Formula 

cct.r 

cci  r, 

car  ccif, 

CCIF  CCIF, 

Molecular  Weight 

137  37 

120  92 

187  38 

170  93 

Bailing  Point  at  1  atm 

*C 

23  82 

-  29  79 

47  57 

3  77 

•r 

74  87 

-21  62 

11763 

38  78 

Proofing  Point 

*c 

-111 

-  158 

-35 

-  94 

*F 

-  168 

-252 

-31 

-  137 

Critical  Temperature 

•c 

138  0 

1120 

214  1 

145  7 

*r 

388  4 

233  6 

4174 

294  3 

Critical  Prtnuro 

atm 

43  5 

40  5 

33  7 

37  2 

fba/aq  In  aba 

639  5 

596  9 

495 

473  2 

Critical  eeiume 

cc /mof 

247 

217 

325 

293 

ce  ft/lb 

0  0289 

00287 

0  0278 

0  0275 

Critical  Oenaity 

|/c« 

05*4 

0  558 

0  576 

0  582 

fba/cu  ft 

34  6 

34  8 

36  0 

36  32 

Density.  liquid 

f'ct 

1  476 

1311 

1  565 

1  456 

at  25“C  (77*F) 

fbl'CU  ft 

92  14 

81  84 

97  69 

90  91 

at  54  S*C  (130T) 

t'et 

1  40 

1  19 

1  49 

1  36 

Iba/ce  ft 

87  40 

74  28 

93  01 

84  90 

Density.  Sai  d  Vapor 

1/1 

5  88 

6  33 

738 

7  83 

at  Soiling  Point 

iba/cu  ft 

0  367 

0  395 

0  461 

0  489 

S pacific  Heat.  liquid 

(H*at  Capacity) 

»!'(«>  CC) 

0  208 

0  232 

0218 

0  243 

at  25*C  (77*  F) 

8tu/(Ib)  (*F) 

Specific  Hoat.  Vapor, 

at  Const  Pressure  (1  aim) 

«•!/(»>  r  o 

0 142  •  38*C 

0,4Z  #  uoo-n 

0,61  •<uo-n 

at  25’C  (77*F> 

•tu/(ib)  rn 

0  145 

0  170 

Specific  Heat  Ratio  of  Vapor 

at  25*C  and  1  atm 

Cp/Cv 

I  137  • 

1 IJ/  •  (100*0 

1.137 

10,0  »,mo-o 

1  0R4 

Hot!  of  Vaporiration 

cil/g 

43  10 

39  47 

35  07 

32  51 

at  Bolling  Point 

Btu/lb 

77.51 

71  04 

63  12 

58  53 

Thermal  Conductivity  at  ZS"C  (77*F) 

Btu -  (hr )  (ft)  <*F) 
liquid* 

0  0506 

0  0405 

0  0434 

0  0372 

Vapor  (1  atm) 

000451 

0.0057S 

o  cn  14 

0  0060 

(Data  from  ASHRAC  In  moat  cam) 

(0  5  atm) 

Viscesity  at  Z5*C  <77*F) 

liquid 

eentlpoue 

0  430 

0214 

0.68 

0  36 

Vapor  (1  atm) 

cantipoiaa 

00107 

00123 

OHIO 

00112 

fOata  from  ASHRAE  In  moit  cates) 

(0  t  atm) 

Surface  Tenjion  at  25*C  (77*F) 

18 

9 

173 

12  > 

dynes/cm 

Retractive  Indot  ol  Liquid 

1  374 

1.287 

1  354 

1  288 

at  Z5*C  ( 7  7  *  F) 

Relative  Dielectric  Sl'ength  of  Vapor  at  1  atm  1 

3  71 

2  46 

3  9  (0  44  atm) 

3  34 

and  75*C  ( 77* F)  (nitrogen 

- 1) 

Dielectric  Conatant 

Liquid 

2  28  #  29*C 

1  0032  c 

2  41  9  25*C 

2  25  »  ?5*C 

Vaoor  (1  atm)‘« 

1  0036  •  24*C9> 

I  004)  9  26  8*C 

lolpMUty  of  "Frron"  »n  Water 

0017 

at  I  atm  and  25'C  (77*F> 

•i% 

Oil 

0  028 

(Sal  n  Prea) 

0013 

Solubility  of  Water  in  "freon" 

at  25*C  <77*0 

■IN 

0  011 

0  009 

0  011 

0  009 

Solubdlty  Parametar  (  ) 

75 

61 

72 

62 

Been  Butanol  Valeo  (88) 

60 

18 

32 

12 

loaicity* 

Threshold  limit 

W"  <>'»l 

1000 

1000 

1000 

1000 

Valeo  (TIV) 

5800 

4950  | 

7600  [ 

7000 
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VII.  RECOMMENDATIONS 


1.  To  prevent  considerable  temperature  variations  along 
and  around  the  boiling  tube,  a  shrink-fit  process 
must  be  applied  in  order  to  have  more  uniform  contact 
pressure  between  the  copper  sleeve  and  the  boiling 
tube . 

2.  To  determine  the  effect  of  contact  resistance 
(between  copper  sleeve  and  boiling  tube)  on  the 
nucleate-boiling  heat - trans f er  performance  of  a  hori¬ 
zontal  tube,  a  series  of  experimental  investigations 
should  be  performed  with  different  contact  pressures. 

3.  To  maintain  constant  electrical  power  input  into  the 
boiling  tube,  a  voltage  regulator  should  be  added  to 
the  input  line. 

4.  To  keep  the  temperature  of  the  water- ethylene  glycol 
mixture  at  a  constant  value,  the  1/2-ton  capacity 
R-12  refrigeration  unit  should  be  replaced  with  a 
2- ton  unit . 

5.  The  positive-displacement  pump  should  be  replaced 
with  a  centrifugal- type  pump,  in  order  to  deliver 
lower  flow  rates  without  reaching  an  electrical  over¬ 
load  limit . 

6.  To  check  the  saturation  temperature  of  the  boiling 
liquid  corresponding  to  the  measured  saturation  pres¬ 
sure,  a  more  accurate  pressure  transducer  should  be 
used. 

7.  The  active  length  of  the  cartridge  heater  should  be 
carefully  measured  or  obtained  from  the  manufacturer. 
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contact  resistance  present  in  the  boiling  tube,  or  to 
the  errors  in  Henrici’s  wall  temperature 
measurement  s . 
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VI.  CONCLUSIONS 


1.  An  experimental  apparatus  to  study  the  nucleate- 
boiling  heat-transfer  performance  of  R-114  from  a 
single  horizontal  tube  has  been  designed,  constructed 
and  instrumented. 

2.  An  analysis  scheme  to  establish  the  nucleate  boiling 
heat- transfer  coefficient  of  R-114  has  been  outlined, 
and  data  acquisition/ reduction  programs  have  been 
written  and  tested. 

3.  The  experimental  apparatus  was  successfully  operated 
in  the  saturation  temperature  range  of  0  °C  to  20  °C 
and  in  the  heat  flux  range  of  1  to  95  kW/m2 . 

4.  The  test  runs  performed  on  different  days  on  the 
apparatus  showed  repeatabilty  within  ±  1  percent, 
revealing  successful  operation  of  the  apparatus. 

5.  The  uncertainty  in  wall  superhet  (about  20%)  drasti¬ 
cally  .affected  the  boiling  heat  -  transfer  coefficient. 
This  can  be  attributed  to  the  large  temperature  vari¬ 
ations  measured  along  and  around  the  boiling  tube 
because  of  a  non-uniform  contact  pressure  between  the 
copper  sleeve  and  the  boiling  tube. 

6.  The  data  taken  on  this  experimental  apparatus  have 
been  compared  with  the  experimental  correlation  for 
refrigerant  and  refrigerant-oil  mixtures  developed  by 
Chongrungreong  and  Sauer  [Ref.  7].  The  present  data 
agreed  to  within  5  percent  with  this  correlation. 

7.  The  current  data  were  also  compared  with  Henrici's 
[Ref.  10]  experimental  nucleate  boiling  data  (from 
smooth  horizontal  copper  tube  in  R-114).  The  present 
data  lie  as  much  as  50  percent  below  Henrici's  data. 
This  disagreement  may  be  attributable  to  the  thermal 
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Figure  5.11  The  Effect  of  Pressure  on  Boiling  Performance. 
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the  Boiling  Performance  with  the 
Chongrungreong  and  Sauer  (Ref.  7) 


Figure  5.7  Reproducibility  of  the  Heat-Transfer 


Figure  5.6  Reproducibility  of  the  Apparatus. 


Photograph  of  the  Boiling  Tube  During  High  Heat  Flux  Opera 


a/A  -  92700  W/S- 
0 /A  =  2000  W/u 


PRESSURE- ENTHALPY  DIAGRAM 


CHANG! 
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APPENDIX  D 

THERMOCOUPLE  CALIBRATION 


A.  EQUIPMENT  USED 

The  equipment  used  in  thermocouple  calibration  is  shown 
in  Figure  D.l,  and  a  brief  description  of  each  component  is 
given  below: 


1.  Thermocouple  Wire 


Type-T  (copper- constantan)  Teflon- coated  wire  of 
0.254  mm  (0.01  in.)  in  diameter  was  used  for  all  thermocou¬ 
ples  . 


2.  Calibration  Bath 


A  Thermos  flask  was  used  as  the  calibration  bath.  In 
order  to  maintain  an  isothermal  temperature  distribution,  a 
motor-driven  mixer  was  used.  Observations  during  the  cali¬ 
bration  procedure  showed  bath  temperature  fluctuations  to  be 
±  0.002  K. 


3 .  Thermocouple  Readout 


A  Hewlett  Packard  3497A  automatic  data  acquisition/ 
control  system  and  a  Hewlett  Packard  9826  computer  were  used 
to  read,  analyze  and  record  the  calibration  data. 


Reference  Temperature 


A  Hewlett  Packard  2804A  quartz  thermometer  was  used 
to  measure  the  bath  temperature.  This  quartz  thermometer  had 
a  resolution  of  0.0001  K,  while  the  manufacturer-guaranteed 
accuracy  was  better  than  ±  0.03  K. 


B.  PREPARATION  FOR  CALIBRATION 


1 .  Thermocouple  Preparation 

A  total  of  five  thermocouples  were  prepared  for  the 
calibration  run.  Two  thermocouples  were  made  from  the 
beginning  of  a  spool  of  copper- constantan  wire;  two  were 
made  from  the  end;  and  one  was  made  at  the  mid  section  using 
the  following  procedure: 

•  The  Teflon  insulation  was  removed  for  a  length  of  about 
4  mm  from  one  end  of  aim  long  piece  of  wire  and  a 
thermocouple  bead  was  made  using  a  Dynatech  Corporation 
thermocouple  welder. 

•  The  other  end  of  the  thermocouple  wire  was  connected  to 
the  data  acquisition  and  control  unit  through  a  junc¬ 
tion  box. 

2 .  Computer  Program 

A  short  computer  program  (TCAL)  was  written  to 
accept  the  thermocouple  readings  through  the  data  acquisi¬ 
tion  system  and  the  bath  temperature  through  the  digital* 
quartz  thermometer.  A  listing  of  the  computer  program,  TCAL 
is  provided  in  Appendix  G.  This  program  prints  all  data  and 
the  discrepancy  (i.e.,  the  Quartz  thermometer  reading  minus 
the  thermocouple  reading)  as  well  as  it  stores  the  data  on  a 
computer  disk.  In  order  to  convert  the  e.m.f.  values  to 
temperature,  the  manufacturer's  conversion  equation  (D.l) 
was  used  in  the  computer  program  TCAL: 


T  =  a0  +  a:El  +  a2E2  +  a3E5  (D.l) 

+  a4E<*  +  a5E5  +  a6Es+a7E7 
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where : 


T  =  temperature  (°C) 
ao=  0.100860910 
ax=  25727.94369 
a2  =  -767345.8295 
a3  =  78025595.81 
a4=  -9247486589 
a5  =  6 . 97688E+ 11 
a6=  -  2 . 66 192E  + 13 
a7=  3 . 94078E+ 14 

E  =  thermocouple  reading  (volts) 


CALIBRATION  PROCEDURE 

Since  the  temperature  measurements  during  this  investi¬ 
gation  ranged  from  25  °C  to  70  °C,  thermocouple  cali¬ 
bration  was  performed  in  this  region.  First,  hot  water 
(at  about  75  °C)  was  added  to  the  thermos  flask  cali¬ 
bration  bath.  Following  proper  mixing,  all  thermocouple 
readings  were  recorded  through  the  data  acquisition 
system  and  the  bath  temperature  was  measured  by  the 
quartz  thermometer.  To  obtain  lower  temperatures, 
small  quantities  of  cold  water  were  gradually  added  to 
the  bath. 

The  discrepancy  (i.e.,  the  Quartz  thermometer  reading 
minus  the  thermocouple  reading)  was  plotted  against  the 
thermocouple  reading  as  shown  in  Figure  D.2.  The 
second-order  polynomial  curve  (D.2)  shown  in  this 
figure  was  generated  using  the  data  from  all  5 
thermocouples . 


~  ZXAS&&  *&&&&  m&&&: 


DCP  =  - 6 . 7422934E- 02  +  9 . 0277043E-03  T  -  (D.2) 

9. 3259917E-05  T2 

where : 

T-  =  thermocouple  reading 

(from  equation  D.l)  (°C) 

DCP  =  discrepancy  (K) 

Thus,  the  corrected  temperature  values  were  obtained, 
using  equation  (D.3): 


Tb  =  DCP  +  T  (D.3) 

where : 

=  actual  calculated  temperature  (°C) 


Schematic  of  the  Calibration  Devices. 


Thermocouple  Calibration  Curve 


APPENDIX  E 

DATA  REDUCTION  PROGRAM 


I  noo  »  PILE  NAME :  DRPR1 
1010!  DATE:  October  19.  1984 

1020!  REVISED:  November  1.  1984 

1090! 

1040  BEEP 

1050  PRINTER  IS  1 

1060  PRINT  USING  "4X .'"'Select  option: . 

1070  PRINT  USING  "6X.'"‘0  Taking  data  or  re-processinq  previous  data . 

1080  PRINT  USING  "6X  . 1  Plotting  data . 

1090  INPUT  Idp 

MOO  IF  ldp-0  THEN  CALL  Main 

1110  IF  Idp*1  THEN  CALL  Plot 

1120  END 

1 1 30  SUB  Main  • 

1140  CUM  /Cc7  C<7> 

1 1 50  DIM  Emf<t2>.T(l2).DIa<2> ,D2a<2> .Dia<2> ,Doa( 2) .La(2) .Lua(2) ,Kcua(2> 
1160  DATA  0.10086091.25727.94369,-767345.8295.78025595.81 
1170  DATA  -9247486589. 6. 97G88E  +  1  1  , -2. 66192EM3. 3. 94078EM4 
1180  READ  C<»> 

I  1 90  FRINTER  IS  701 
1200  CLEAR  709 
1210  BEEP 

1220  INPUT  "ENTER  MONTH.  DATE  AND  TIME  <MM : DD :HH :MM : SS)" , DateS 
1230  OUTPUT  709 : "TO" : DateS 
1240  OUTPUT  709:"TD" 

1250  ENTER  709:Date$ 

1260  PRINT 
1270  PRINT 
1280  PRINT 

1290  PRINT  "  Month,  date  and  time  :":Date$ 

1300  PRINT 

1310  PRINT  USING  "1  OX , ""NOTE :  Program  name  :  DRP . 

1320  BEEP 

1330  INPUT  "ENTER  DISK  NUMBER". Dn 

1340  PRINT  USING  "16X,""Disk  number  ■  "".ZZ":Dn 

1 350  PEEP 

1360  INPUT  "ENTER  INPUT  MODE  <0-3054A, 1 -FILE)"  .  Iro 
1370  IF  Im-0  THEM 
I  380  BEEP 

1390  INPU1  "GIVE  A  NAME  F1JR  THE  RAW  DATA  FIIE".D2  Tile? 

1400  PRINT  USING  " 1 6X . ""New  file  name:  , 1 4A" ; D2_fileS 

1410  CREATE  BOAT  D2_file$,20 
1420  ASSIGN  »FHe2  TO  D2_file$ 

I '130! 

1440!  DUMMY  FILE  UNTIL  Nrun  KNOWN 
1490  D1  f i IeS-"DUMHY" 

1460  CREAIE  BDAT  D1  fl leS,20 
1470  ASSIGN  (SFllel  To  D1_fileS 
1480  0UIP1JT  SFllel  :Date$ 

1490  BEEP 

1500  INPUT  "GIVE  A  NAME  FOR  THE  PLOT  FIIE".P  fileS 
1510  CREATE  BDAT  P_f t leS.S 
1520  ASSIGN  jiPlot  TO  P.fileS 
1530  BEEP 

1540  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  ICS  < 0-DEFAULT >", Id tc 
1550  IF  Idtc'O  THEN 
1560  Ldtcl-0 
1570  Ldtc2-0 

1580  PRINT  USING  "16X,""No  defective  TCs  exist""" 

1530  END  IF 


1600  IF  Idtc-t  THEN 
1610  BEEP 

1620  INPUT  "ENTER  DEFECTIVE  TC  LOCATION" . Ldtc 1 

1630  PRINT  USING  "16X,""TC  is  defective  at  location  "".D":Ldtc1 

1 640  Ldtc2-0 

1650  END  IF 

1660  IF  Idtc-2  THEN 

1670  BEEP 

1680  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS" .Ldtcl .Ldtc2 

1690  PRINT  USING  "16X,""TC  are  defective  at  locations  B  .  4X  ,  D"  :  Ldtcl  .  Ldtc2 

1700  END  IF 

1710  IF  Idtc>2  THEN 

1720  BEEP 

1730  PRINTER  IS  1 

1740  BEEP 

1750  PRINT  "INVALID  ENTRY" 

1760  PRINTER  IS  701 
1770  GOTO  1530 
1780  END  F 

1790  OUTPUT  ?File1 :Ldtc1 ,Ldtc2 
1800!  Im-1  option 
1810  ELSE 
1820  BEEP 

1830  INPUT  "GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" ,D2_f i le$ 

1840  PRINT  USING  "16X.""01d  file  name:  "" , 1 4A" : D2_f i leS 
1850  ASSIGN  aFile2  TO  D2_file$ 

1860  ENTER  9File2:Nrun 
1870  ENTER  ?File2;Dold$ 

1 8P0  BEEP 

1890  INPUT  "GIVE  A  NAME  FOR  PLOT  F ILE" ,P_f i le$ 

1900  CREATE  BDAT  P_fileS.5 
1910  ASSIGN  SPlot  TO  P_fileS 

19 20  PRINT  USING  "16X,""This  data  set  taken  on  :  . 1 4A” : Do ld$ 

1930  ENTER  »File2:Ldtc1 .Ldtc2 

1940  PRINT  USING  " 1 6X . ""Thermocoup les  were  defective  at  locat ions : "" . 2< 3D . AX) " ; 

Ldtcl , Ld  tc2 

1 950  ENTER  9File2;Itt 

1960  END  IF 

1970  PRINTER  IS  I 

1980  IF  Im-0  THEN  .' 

'990  BEEP 

2000  PRINT  USING  "AX. ""Select  tube  type . 

2010  PRINT  USING  "6X , ""0-Smooth  4  inch  Ref . 

2020  PRINT  USING  "6X  ."" 1 -Smooth  4  inch  soft  Cu . 

2030  PRINT  USING  "6X . ""2-Smooth  8  inch  soft  Cu”"" 

2040  INPUT  Itt 

2050  IF  1 1 1 > 2  THEN 
2060  BEEP 

2070  PRINT  "INVALID  ENTRY" 

2080  GOTO  2000 

2090  END  IF 

2100  OUTPUT  AFi lei : 1 1 1 

21 10  END  IF 

2120  PRINTER  IS  701 

2130  PRINT  USING  "16X.""Tube  Type  is  :  "".B":Itt 
2140  BEEP 

2  150  INPUT  "ENTER  OUTPUT  VERSION  ( 0-LONG . 1 -SHORT  > "  ,  Iov 
2160! 

2170!  Dl-Diameter  at  thermocouple  positions 
2180  DATA  .01 1 125. .01 1 1 125. .01 143 
2190  READ  Dla<*> 
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2200  DI-Dla(Itt) 

2210! 

2220!  D2-Diameter  of  test  section  to  the  base  of  fins 
2230  DATA  . 01 5875 .. 01 5875 ., 0 1 5875 
2240  READ  D2a<-> 

2250  D2*D2a< 1 1 1  > 

2260! 

2270!  Di-Inside  diameter  of  unenhanced  ends 
2280  DATA  . 0 1 27 , . 01 27 . . 0 1 27 
2290  READ  Dia<*> 

2300  Di»Dia< Itt) 

2310! 

2320!  Do=0utside  diameter  of  unenhanced  ends 
2330  DATA  .  0 1 5875 .. 01 5875 0 1 5875 
2340  READ  Doa<»> 

2350  Do-Doa<Itt> 

2360! 

2370!  L-Length  of  enhanced  surface 
2380  DATA  . 1016, . 1016. .2032 
2390  READ  La<-> 

2400  L=La< 1 1 1  > 

2410! 

2420!  Lu*Length  of  unenhanced  surface  at  the  ends 
2430  DATA  .  0254 0254 .. 0762 
2440  READ  Lua< »  > 

2450  Lu*Lua<Itt> 

2460! 

2470!  Kcu-Thermal  Conductivity  of  tube 
2430  DATA  398.344,344 
2490  READ  Kcua(-) 

2500  Kcu-Kcua< I t t ) 

2510  A-PI*<Do  2-Dl'2)/4 

2520  P-PI*Do 

2530  J-l 

2540  Sx-0 

2550  Sy-0 

2560  Sxs-0 

2570  Sxy-Q 

2580  Repeat :  ! 

2590  IF  Im-0  THEN 

2600  ON  KEY  0.15  RECOVER  2530 

2610  PRINTER  IS  I 

2620  PRINT  USING  ”4X  .  ""SELECT  OPTION . 

2630  PRINT  USING  "9X  ,  ""0  =  TAKE  DATA . 

2640  PRINT  USING  "9X,""1=SET  HEAT  FLUX . 

2650  PRINT  USING  "9X , ”"2=SE T  Tsat . 

2660  PRINT  USING  "4X  . ""NOTE :  KEY  0  >  ESCAPE . 

2670  BEEP 

2680  INPUT  Ido 

2690  IF  Ido-0  THEN  3630 

2700! 

27 !  0 !  LOOP  TO  SET  HEAT  FLUX 

2720  IF  Ido-1  THEN 

2730  OUTPUT  709:”AR  AF62  AL63  VR5" 

2740  BEEP 

2750  INPUT  "ENTER  DESIRED  0dp".Dqdp 

2760  PRINT  USING  "4X , ""DESIRED  Qdp  ACTUAL  Qdp . 

2?70  Err-1000 

2780  FOR  1-1  TO  2 

2790  OUTPUT  709:"AS  SA" 

2800  Sum-0 
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2b  10  T  OR  J 1  -  I  ro  5 
2820  ENTER  709:E 
2830  Sum-Sum+E 
2840  NEXT  Ji 

3850  IF  1-1  THEN  Vo  It -Sum-5 
2860  IF  1-2  THEN  Amp-E 
2870  NEXT  I 

2880  Aqdp-Vol t*Amp/ (PI-D2-L) 

2890  IF  ABS< Aqdp-Dqdp ) >Err  THEN 

2900  IF  Aqdp>3qdp  THEN 

2910  BEEP  4000,  .2 

2320  BEEP  4000.. 2 

2930  BEEP  4000. .2 

2940  ELSE 

2950  BEEP  250.  .2 

2960  BEEP  250.  .2 

2970  BEEP  250.. 2 

2980  END  IF 

2990  PRINT  USING  "4X . MZ . 3DE  , 2X , HZ . 3DE" : Dqdp . Aqdp 

3000  WAIT  2 

3010  GOTO  2780 

3020  ELSE 

3030  BEEP 

3040  PRINT  USING  "4X  ,MZ  .  3DE  ,2X .HZ . 3DE" : Dqdp . Aqdp 

3050  Err-500 

3060  WAIT  2 

3070  GOTO  2780 

3080  END  IF 

3090  END  IF 

3100! 

3!  10!  LOOP  TO  SET  Tsat 

3120  IF-  Ido-2  THEN 

3130  OUTPUT  709:"AR  AF33  AL33  VR5" 

3140  OUTPUT  709:  "AS  SA" 

3150  BEEP 

3160  INPUT  "ENTER  DESIRED  Tsaf'.Dtld 

3170  PRINT  USING  "4X.""D  Tsat  A  Tsat  DiFf  Psat 

3180  01d-0 

3190  ENTER  709 : E 1  lq 

3200  At ld-FNTvsv<  E 1 iq  > 

3210  Psat-FNPsat  <  At  Id) 

322 0  IF  ABS<Atid-0tld>>.2  THEN 

3230  IF  AtlrODtld  THEN 

3240  BEEP  4i.  30.  .2 

3250  BEEP  4000.  .2 

3260  BEEP  4000.  .2 

3270  ELSE 

3280  BEEP  250. .2 

3290  BEEP  250.  .2 

3300  BEEP  250.. 2 

3310  END  IF 

3320  Err-Atld-Old 

3330  Old-At Id 

3340  PRINT  USING  "4X . 4< HDD . DD . 3X X" ; 0 t 1 d . A t Id . E r r , Psat 
3350  WAIT  2 
3360  GOTO  3190 
3370  ELSE 

3380  IF  ABS< At Id-D t Id) > . 1  THEN 
3390  IF  At  ld>Dt Id  THEN 
3400  BEEP  3000.. 2 
3410  BEEP  3000. .2 
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3420  ELSE 
3430  BEEP  800.. 2 
3440  BEEP  300.  .2 
3450  END  IF 
34G0  E  r  r -A  t Id-Old 
3470  Old-Atld 

3480  PRINT  USING  "4X  .  4(MDD  .  DD  .  3X )  " ;  D  t  Id  .  At  1  d  ,  Er  r  .  Psat 

3490  WAIT  2 

3500  GOTO  3190 

3510  ELSE 

3520  BEEP 

3530  Err-Atld-Old 

3540  Old-Atld 

3550  PRINT  USING  "4X . 4 (MDD . DD . 3X >" ; Dt Id . At  Id .Er r . Psat 

3560  WAIT  2 

3570  GOTO  3190 

3580  END  IF 

3590  END  IF 

3600  END  IF 

3610!  ERROR  TRAP  FOR  Ido  OUT  OF  BOUNDS 

3620  IF  Ido>2  THEN 

3630  BEEP 

3640  GOTO  2620 

3650  END  IF 

3660! 

3670!  TAKE  DATA  IF  Im-0  LOOP 
9680  BEEP 

3690  INPUT  "ENTER  BULK  OIL  %".Bop 
3700  OUTPUT  709 : "AR  AF25  AL36  VR5" 

3710  FOR  1-1  TO  12 
3720  OUTPUT  709:"AS  SA" 

3730  Sum-0 
3740  FOR  Ji-1  TO  20 
3750  ENTER  709:E 
3760  Sum-Sum+E 
3770  NEXT  Jt 
3780  Emf  ( I>-Siun/20 
3790  NEXT  I 

3800  OUTPUT  709 : "AR  AF62  AL63  VR5" 

3810  FOR  1-1  TO  2 
3820  OUTPUT  709: "AS  SA" 

3830  Sum-0 
3840  FOR  Ji-1  TO  20 
3850  ENTER  709  :E 
3860  Sum=Sum+E 
3870  NEXT  Ji 

3880  IF  1-1  THEN  Vr-Sum/20 
3890  IF  1-2  THEN  Ir-Sum/20 
3900  NEXT  I 
3910  ELSE 

3920  ENTER  «F i le2 :Bop . To IdS ,Emf < * >  . Vr . I r 

3930  END  IF 

3940! 

3950!  CONVERT  emf'S  TO  TEMP .VOLT .CURRENT 

3960  Tua-0 

3970  FOR  I- I  TO  12 

3980  IF  Idtc>D  THEM 

3990  IF  I  -Ld  tc  I  OR  I-Ldtc2  THEN 

4000  T<  I  > - -99 . 99 

4010  GOTO  4060 

4020  END  IF 
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4040  T  ( I  )*FNT>/sv(Emf  <  I  >  > 

4050  IF  I  <9  THEN  Tua* Tua+T < I > 

4060  NEXT  I 

4070  Twa*Tij.ia/<8-Idtc) 

40F0  Tld*T<9> 

4000  Tv*-  <  T  <  1  0  )  +  T  <  1 1  ) )  /2 
4  1 00  Isump*  T < 12) 

4110  Amp  *  I r 
4120  Vo  1 t*Vr*25 
4130  Q  =  Vo 1 t*  Amp 
4140  IF  Itt-0  THEN 
4150  Kcu*FNKcu< Tua) 

4160  ELSE 

4170  Kcu*Kcua< I t t > 

4130  END  IF 
4190! 

4200!  FOURIER  CONDUCT  ION  EQUATION  HUH  CONTACT  RESISTANCE  NEGLECTED 
4210  Tu*Twa-G*L0G<D2/D1  >/<2*PI*Kcu»L> 

4220  The  tab* Tu- T 1 d 
4230  IF  Ihe tab<0  THEN 
4240  BEEP 

4250  INPUT  "TWALL<T$AT  <0*CONTINUE.  1-END) ".lev 

4260  IF  lev*0  THEN  GOTO  2590 

4270  IF  lev* 1  THEN  5030 

4280  END  IF 

4290! 

4300!  COMPUTE  VARIOUS  PROPERTIES 
4310  Tfilm*FNTf ilm<Tw,Tld) 

4320  Rho*FMRho<Tf lira) 

4330  Mu*FNMu(Tf tlm) 

4  340  K-FNK <  T f i lm ) 

4350  Cp”FNCp<  T f i lm) 

4360  Be ta*FN8e tat Tf i lm) 

4370  N i *Mu./Rho 
4380  A  lpha  =  K/ ( Rho*Cp ) 

4390  Pr  *N i /Alpha 
4400  Psat*FNPsat<Tld) 

4410! 

4420!  COMPUTE  HEAT  TRANSFER  COEFFICIENT 
4430  Hbar  *  1 90 

4440  Fe= < Hbar*P/ < Kcu*A ) )  . 5*Lu 
4450  Tanh=FN Tanh < Fe ) 

4a60  The ta* The  tab* T anh/Fe 
4470!  PRINTER  IS  701 

4490!  PRINT  USING  "4X . 7 < 1 X . MZ . 3DE > " : Hba r . Fe . Tanh . The  tab . The ta . Be ta. N i 
4490  X» * < 9 . 8 1 “Be ta* The tab»Do  3* Tanh/ ( Fe*N i * A  1 pha> )  .166667 
4500  Yy  =  <  1  M  ,559/F'r)  <9/16>>  <  8/27  > 

4510  Hbarc-K/Do*( .6+ . 387*Xx/Yy )  2 

4520  IF  ABS < < Hba r - Hbarc ) /Hbarc ) > . 00 1  THEN 

4530  Hbar *( Hbar “Hbarc >*. 5 

4540  GOTO  4440 

4550  END  IF 

4660  Q 1  * <Hbar *P»Kcu*A )  . 5* The  tab* Tanh 

4570  Qc*Q-2*Q  1 

4580  As  =  P  I  *D2*L 

4590  Qdp*Qc/A3 

4600  H t ube*Qdp/ T he t ab 

4610! 

4620!  RECORD  TIME  OF  DATA  TAKING 
4630  IF  lm*0  THEN 
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4640  OUTPUT  709;"TD" 

4650  ENTER  709:ToldS 

4660  END  IF 

4670! 

4680!  OUTPUT  DATA  TO  PRINTER 
4690  PRINTER  IS  701 
4700  PRINT 
4710  IF  IowO  THEN 

4720  PRINT  USING  "1 IX  .""Data  Set  Number  -  "". DDD . 2X . "“Bulk  Oil  '/.  *  ,DD .  D  .5X .  1 
4A” : J.Bop . Told? 

4730  PRINT 

4740  PRINT  USING  "11X.""TC  No:  1  2  3  4  5  6  7 

6 . 

4750  PRINT  USING  " 1 ! X . ""Temp  : "" ,8< 1 X .HDD . DD> " : T < 1 ) . T<2> . T<3> . T<4> , T < 5 > . T<6> . T< 
7  > . T<8) 

4760  PRINT  USING  "11X.""  Twa  Tliqd  Thetab  Tvapr  Psat  Tsump . 

4770  PRINT  USING  "1 1 X ,2<MDD . DD . 1 X) .HZ . DDDE . IX ,2< 1 X .MDD.DD) ,2X .MDD .  D" ; Twa. Tld . Th 
e  tab , Tv . Psat . Tsump 

4780  PRINT  USING  “1 IX,""  Htube  Odp . 

4790  PRINT  USING  "1 1 X ,MZ ,3DE , 1 X ,MZ . 3DE" : Htube .Qdp 

4800  END  IF 

4810  IF  Iov-1  THEN 

4820  IF  J-l  THEN 

4830  PRINT  USING  "11X.""RUN  No  Oil*/.  Tsat  Htube  Qdp  Thetab . 

4840  END  IF 

4850  PRINT  USING  " 12X . DDD . 4X , DD ,2X .MDD .DD . 3<1 X ,MZ . DDDE)" : J . Bop , T Id .Htube . Qdp . Th 
etab 

4860  END  IF 
4870  IF  Im-0  THEM 
4880  BEEP 

4890  INPUT  "OK  TO  STORE  THIS  DATA  SET  <  I  -  Y  .  0-N>  ?"  .  Olt 
4900  END  IF 

4910  IF  Ok  - 1  OR  Im-I  THEN  J-J+l 

4920  IF  Ok  - 1  AND  Im-0  THEN  OUTPUT  AFi lei : Bop . To IdS . Emf < » > . Vr . Ir 
4930  IF  Im-1  OR  Ok  - 1  THEN  OUTPUT  APlot :Qdp . Thetab 
4940  IF  Im-0  THEN 
4950  BEEP 

4960  INPUT  "WILL  THERE  BE  ANOTHER  RUN  < 1 -T . 0-N>?" .Go_on 
4970  Nrun-J 

4980  IF  Go_on<> 1  THEN  5030 
4990  IF  Go_on- 1  THEN  Repeat 
5000  ELSE 

5010  IF  J<Nrun+l  THEN  Repeat 

5020  END  IF 

5030  IF  Im-0  THEN 

5040  BEEP 

5050  PRINT 

5060  PRINT  USING  "10X.""N0TE:  "".ZZ.""  data  runs  were  stored  in  file  "".I0A":J- 
1 .D2_f i le$ 

5070  ASSIGN  siFilel  TO  « 

5080  OUTPUT  @File2:Nrun-1 
5090  ASSIGN  «Ftlel  TO  D1  file? 

5100  ENTER  AF l  lei : Date* .Ldtcl ,Ldtc2 . 1 1 1 
5 M 0  OUTPUT  9File2:Date*. Ldtcl ,Ldtc2.Itt 
5120  FOR  1-1  TO  Nrun-1 
5130  ENTER  AF  i  le  1  ; Bop .  folds ,Emf < • > , Vr . I r 
5140  OUTPUT  AF  i  le2:  Bop  .  Told*.Emf<»)  .Ur.Ir 
5150  NEXT  I 

5160  ASSIGN  AF i  lei  TO  - 
5170  PURGE  "DUMMY" 

5180  END  IF 


.ZZ.""  data  runs  were  stored  in  file  "".I0A":J- 


ji3u  Dtcr' 
t>200  PR1NI 

5210  PRINT  USING  "1  OX. ""NOTE:  "".ZZ.""  X-Y  pairs  ye  re  stored  in  plot  data  file 
"" , 1 0ft" : J- 1 ,P_f ile$ 

5220  ASSIGN  9F!le2  TO  - 
5230  ASSIGN  SPlot  TO  • 

5240  BEEP 

5250  INPUT  "LIKE  TO  PLOT  DATA  <  1  *Y  .0-N>r\0k 

5250  IF  0k  - 1  THEN 

5270  CALL  Plot 

5280  END  IF 

5290  SUBEND 

5300! 

5310!  CURVE  FITS  OF  PROPERTY  FUNCTIONS 
5320  DEF  FNKcu(T) 

5330!  OFHC  COPPER  250  TO  300  K 

5340  Tk-T+273.15  !C  TO  K 

5350  K-434- . I  1 2-Tk 

5360  RETURN  K 

5370  FNEND 

5380  DEF  FNMu(T) 

5390!  170  TO  360  K  CURVE  FIT  OF  VISCOUSITY 

5400  Tk-T+273.15  !C  TO  K 

5410  Mu“EXP< -4 . 4636+ < 1011 . 47/Tk ) )* 1 . OE-3 

5420  RETURN  Mu 

5430  FNEND 

5440  DEF  FNCp<T> 

5450!  180  TO  400  K  CURVE  FIT  OF  Cp 

5460  Tk-T+273.15  !C  TO  K 

5470  Cp-. 40188+1 .65007E-3»Tk +1 .51494E-6«Tk  2-6 .67853E- 1 0«Tk  3 

5480  RETURN  Cp 

5490  FNEND 

5500  DEF  FNRho<  T> 

5510  Tk -T+273 . 15  !C  TO  K 

5520  X- 1  - < 1 .8*Tk/753.95)  !K  TO  R 

5530  Ro-36. 32+61 . 146414-X  < 1/3)+16.418015-X+l 7.476838-X  .5+1 . 1 19828-X'? 

5540  RETURN  Ro 
5550  FNEND 
5560  DEF  FNPr(T) 

5570  Pr-FNCp<T)-FNMu(T)/FNK(T) 

5580  RETURN  Pr 
5590  FNEND 
5600  DEF  FNK<  T) 

5610!  T<360  K  WITH  T  IN  C 
5620  K* . 071  - . 000261 *T 
5630  RETURN  K 
5640  FNEND 
5650  DEF  FNTanM  X) 

5660  P-EXP(X) 

5670  Q-EXP(-X) 

5680  Tanh-(P-0)/(P+0> 

5690  RETURN  Tanh 

5700  FNEND 

5710  DEF  FNTvsv(V) 

5720  COM  7Cc/  C(7) 

5730  T-C<0) 

5740  FOR  1-1  TO  7 
5750  T - T+C < I ) -V  I 
5760  NEXT  I 

5770  T-T-6. 7422934E-2+T*<9. 027 7043E-3-T-( -9. 325991 7E-5>> 
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5780  RETURN  T 

5790  FNEND 

5800  DEF  FNBeta(T) 

5810  Rop-FMRho< T+ . I ) 

5820  Rom-FNRho(T-.l) 

5830  Beta--2/(Rop+Rom)«(Rop-Rom>/ .2 
5840  RETURN  Beta 
5850  FNEND 

5860  DEF  FNTfilm(Tu.Tld> 

5870  T f i  lm-<  Tw  +  Tld)/2 
5880  RETURN  Tfilm 
5890  FNEND 
5900  DEF  FNPsat<  Tc) 

5910!  0  TO  80  deg  F  CURVE  FIT  OF  Psat 
5920  Tf - 1 . 8*Tc+32 

5930  Pa*5 . 945525 + T f  »( . I5352082+Tf «( 1 .4840963E-3+Tf-9.6150671E-6> ) 
5940  Pg-Pa-14.7 

5950  IF  Pg>0  THEN  !  +-PSIG.  — In  Hg 

5960  Psat-Pg 
5970  ELSE 

5980  Psat-Pg-29 . 92/ 14.7 
5990  END  IF 
6000  RETURN  Psat 
6010  FNEND 
6020  SUB  Plot 
6030  DIM  C<9) 

6040  INTEGER  It 
6050  PRINTER  IS  1 
6060  Idv-1 
6070  BEEP 

6080  INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  <  I -Y  ,0-N)?”.Idv 
6090  BEEP 

6100  PRINT  US I MG  "4X. ""Select  Option: . 

6110  PRINT  USING  "4X.""0  q  versus  delta-T . 

6120  PRINT  USING  "4X.""1  h  versus  delta-T . 

6130  PRINT  USING  "4X.""2  h  versus  q . 

6140  PRINT  USING  "4X,""3  h-ratio  versus  delta-T . 

6150  INPUT  Opo 

6160  IF  Opo-3  THEN 
6170  BEEP 

6130  INPUT  "SELECT  TUBE  DIAMETER  (0-. 75, 1-1.0  IN)".Itd 

6190  END  IF 

6200  PRINTER  IS  705 

6210  IF  IdvOI  THEN 

6220  BEEP 

6230  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS". Cx 
6240  BEEP 

6250  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS". Cy 
6260  BEEP 

6270  INPUT  "ENTER  MIN  X-VALUE  (MULTIPLE  OF  10)".Xmtn 
6280  BEEP 

6290  INPUT  "ENTER  MIN  Y-VALUE  (MULTIPLE  OF  IO>".Ymln 

6300  ELSE 

6310  Cy-3 

6320  IF  Opo-O  THEN 

6330  Cx-2 

6340  Xmin*1 

6350  Ymln-1000 

6360  END  IF 

6370  IF  Opo* 1  THEN 

6380  Cx-2 


oodU  Amin-I 

6400  Ymin*100 
6410  END  IF 

6420  IF  Opo-2  THEN 
6430  Cx-3 

6440  Xnin-1000 

6450  Ymin-100 
6460  END  IF 

6470  IF  Opo-3  THEN 
6480  Cx-2 

6490  Cy *3 

6500  Xmin-1 

6510  Ymin*l 

6520  END  IF 

6530  END  IF 

6540  BEEP 

6550  PRINT  "IN:SP1 :  IP  2300.1800.8300.6800:" 
6560  PRINT  "SC  0 . 1 00 , 0 . 1 00 : TL  2.0:" 

6570  Sfx-100/Cx 
6580  SFy*100/Cy 
6590  PRINT  "PU  0.0  PD" 

6600  Nn-9 

6610  FOR  I-l  TO  Cx*l 
6620  Xat-Xmin-10  <I-1> 

6630  IF  I-Cx ♦  I  THEN  Nn-1 

6640  FOR  J-1  TO  Nn 

6650  IF  J-1  THEN  PRINT  "TL  2  0" 

6660  IF  J-2  THEN  PRINT  "TL  1  0" 

6670  Xa-Xat*J 

6680  X-LGT<Xa/Xmin>«Sfx 

6690  PRINT  "PA" :X ,"  ,0  :  XT;" 

6700  NEXT  J 

6710  NEXT  I 

6720  PRINT  "PA  100.0;PUi" 

6730  PRINT  "PU  PA  0.0  PD" 

6740  Nn-9 

6750  FOR  1-1  TO  Cy+1 
6760  Yat-Ymin-10  <  I- 1  > 

6770  IF  I-Cy+1  THEN  Nn-I 

6780  FOR  J-1  TO  Nn 

6790  IF  J-1  THEN  PRINT  "TL  2  0" 

6800  IF  J-2  THEN  PRINT  "TL  1  0" 

6810  Ya-Yat*J 

6820  Y-LRT(Ya/Ymin)-Sfy 

6830  PRINT  "PA  0,";Y."YT" 

6840  NEXT  J 

6850  NEXT  I 

6860  PRINT  "PA  0.100  TL  0  2" 

6870  Nn-9 

6880  FOR  1-1  TO  Cx+1 
6890  Xat-Xmin-10  <I-1> 

6900  IF  I-Cx 1  THEN  Nn-I 

6910  FOR  J-1  TO  Nn 

6920  IF  J-1  THEN  PRINT  "TL  0  2" 

6930  IF  J>1  THEN  PRINT  "TL  0  1" 

6940  Xa-Xat-J 

6950  X-LGT(Xa/Xmin)*Sfx 

6360  PRINT  "PA";X,". 100:  XT" 

6970  NEXT  J 

6980  NEXT  I 


6990  PRINT  "PA  100,100  PU  PA  100.0  PDM 
7000  Nn-9 

7010  FOR  1*1  TO  Cy  +  1 
7020  Yat-Ymtn*10  < I- 1 > 

7030  IF  I-Cy+1  THEN  Nn- I 

7040  FOR  J-1  TO  Nn 

7050  IF  J-1  THEN  PRINT  MTL  0  2" 

7060  IF  J>1  THEN  PRINT  "TL  0  1“ 

7070  Ya-Yat-J 

7080  Y-LGT(Ya/Ymin)»STy 

7090  PRINT  "PD  PA  1 00  , "  .  Y  .  ”YT" 

7100  NEXT  J 

7110  NEXT  I 

7120  PRINT  "PA  100.100  PU” 

7130  PRINT  "PA  0.-2  SR  1  .5.2" 

7140  Ii-LGT<Xmin) 

7150  FOR  1-1  TO  Cx  +  1 
7160  Xa*Xmin« 1 0  ' < I- 1 ) 

7170  X-LGT  <Xa/Xmin)»Sf  x 
7180  PRINT  "PA":X,".0i" 

7190  IF  Ii>-0  THEN  PRINT  "CP  -2 . -2 :LB1 0 ; PR  -2.2:LB";Ii 

7200  IF  I i <0  THEN  PRINT  "CP  -2 , -2 ;LB1 0 : PR  0 .2  :LB" :  Ii 

7210  Ii-Ii+t 

7220  NEXT  I 

7230  PRINT  "PU  PA  0.0" 

7240  Ii-LGT(Ymin) 

7250  Y10-10 

7260  FOR  1-1  TO  Cy+1 

7270  Ya-Ymin-10  (1-1) 

7280  Y-LGT<  Ya/Ymin>*Sf  y 
7290  PRINT  "PA  0,":Y,"" 

7300  PRINT  "CP  -4.-.25;LB10;PR  -2 .2 ;LB" : Ii :"" 

7310  Ii-Ii+1 
7320  NEXT  I 
7330  IF  IdvOI  THEN 
7340  BEEP 

7350  INPUT  "ENTER  X-LABEL" .XlabelS 
7360  BEEP 

7370  INPUT  "ENTER  Y-LABEL" . Y label* 

7380  ELSE 

7390  IF  Opo-O  THEN 

7400  Xlabel$*"Tw-Tsat  <K>" 

7410  YlabelS-"q  <W/m  2>" 

7420  END  IF 

7430  IF  Opo-1  THEN 

7440  X labels-" Tu- Tsat  OO" 

7450  YlabelS-"h  <W/m'2.K>" 

7460  END  IF 

7470  IF  Opo-2  THEN 

7480  XlabelS-"q  (N/m"2)" 

7490  YlabelS-"h  (W/m  2.K)" 

7500  END  IF 

7510  IF  Opo-3  THEM 

7520  XlabelS-"Tw-Tsat  <K)" 

7530  Y labe lS-"h<enh) /h< smooth)" 

7540  END  IF 
7550  END  IF 

7560  PRINT  "SR  1 . 5 . 2 : PU  PA  50.-16  CP" : -LENT X labe 1 S > /2 : "0 : 1.B" : X labe 1 S 
7570  PRINT  "PA  -14.50  CP  0 . " : -LEN< Y labe l$>/2*5/6 : "DI  0 . 1 ;LB" ; Y label® 
7580  PRINT  "CP  0,0  DI" 

7590  Repeat : ! 


/bUU  bfctK 

76(0  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  < 1  - Y  . 0-N>?" .Ok 
7620  IF  Ok  *1  THEN 
7630  BEEP 

7640  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" . D_f i le$ 

7650  ASSIGN  9File  TO  D_file$ 

7660  BEEP 
7670  BEEP 

7680  INPUT  "ENTER  THE  BEGINNING  RUN  NUMBER". Md 
7690  BEEP 

7700  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 

7710  BEEP 

7720  PRINTER  IS  I 

7730  PRINT  USING  "4X. ""Select  a  symbol: . 

7740  PRINT  USING  "4X.""1  Star  2  Plus  sign 

7750  PRINT  USING  "4X.""3  Circle  4  Square 

7760  PRINT  USING  "4X.""5  Rombus . 

7770  PRINT  USING  "4X.""6  R iqht-side-up  triangle . 

7780  PRINT  USING  "4X,""7  Up-side-down  triangle . 

7790  INPUT  Sym 
7800  PRINTER  IS  705 
7810  PRINT  "PU  DI" 

7820  IF  Sym-1  THEN  PRINT  "SM*" 

7830  IF  Sym-2  THEN  PRINT  "SM+" 

7840  IF  Sym-3  THEN  PRINT  "SMo" 

7850  IF  Md> 1  THEN 
7860  FOR  1-1  TO  <Md-1> 

7870  ENTER  9File:Ya.Xa 

7880  NEXT  I 

7890  END  IF 

7900  FOR  1-1  TO  Npairs 

7910  ENTER  9F i le : Ya, Xa 

7920  IF  Opo-1  THEN  Ya-Ya/Xa 

7930  IF  Opo-2  THEN 

7940  0- Ya 

7950  Ya-Ya/Xa 

7960  Xa-Q 

7970  END  IF 

7980  IF  Opo-3  THEN  Ya-Ya/FNHsmooth<Xa. Itd> 

7990  X-LGT  <Xa/Xmin)-Sf  x 
8000  Y-LGT <Ya/Ymln)«Sfy 
8010  IF  Sym>3  THEN  PRINT  "SM" 

8020  IF  Sym<4  THEN  PRINT  "SR  1.4, 2. 4" 

8030  PRINT  "PA" . X . Y . "" 

8040  IF  Sym>3  THEM  PRINT  "SR  1.2. 1.5" 

8050  IF  Sym-4  THEN  PRINT  "UC2 . 4 . 99 . 0 . -8 . -4 . 0 , 0 . 8 , 4 . 0  :  ” 
8060  IF  Sym-5  THEN  PRINT  "UC3. 0 .99 . -3 . -6 , -3.6 .3 .6 .3  . -6 :" 
8070  IF  Sym-6  THEN  PRINT  "UCO .5 . 3 .99.3 . -8 . -6 . 0 .3 .8 : " 

8080  IF  Sym-7  THEN  PRINT  "UCO , -5 .3 .99. -3 ,8 .6 . 0 , -3. -8:" 
8090  NEXT  I 
8100  BEEP 

8110  ASSIGN  9F i le  TO  * 

8120  GOTO  7600 
3130  END  IF 
8140  PRINT  "PU  SM" 

3150  BEEP 

3160  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  ( 1  * Y . 0-N>?" ,Go_on 
8170  IF  Go_on* 1  THEN 
3180  BEEP 

8190  INPUT  "ENTER  LOWER  AND  UPPER  X-L IMITS" . X 1 1  . X lu 


8200  FOR  Xx«0  TO  Cx  STEP  Cx/200 
8210  Xa-Xmin-10  Xx 
8220  IF  Xa<Xl 1  DR  Xa>Xlu  THEN  8290 
3230  Ya-FNPoly<Xa) 

8240  Y-LGT  <Ya/Ymin>*Sfy 
8250  X-LGT<  Xa/Xmln)«Sf  x 
8280  IF  Y<0  THEN  Y-0 
8270  IF  Y> 1 00  THEN  GOTO  8290 
8280  PRINT  "PA" ,X , Y ,"PDM 
8290  NEXT  Xx 
8300  END  IF 

8310  PRINT  "PU  PA  0.0  SP0“ 

8320  SUBEND 

3330  DEF  FNHsmooth< X , I td) 

8340  Hs-FNPoly<X)/X 

8350  IF  I td- 1  THEN  Hs-Hs*. 83347 

8360  RETURN  Hs 

8370  FNEND 

8380  DEF  FNPoly <X> 

8390  Poly--4. 41237 18E+2-X»<6. 81 2391 7E+2-X-3. 74 1S8S3E+2) 

8400  RETURN  Poly 
8410  FNEND 

9420  DEF  FNPvst< Tsteam) 

8430  DIM  K<8> 

8440  DATA  -7.691234564,-26.08023696.-168. 1706546.64.23285504.-1 18. 964622S 
8450  DATA  4 . 1671 1 732 .20 . 9750676 . 1 .E9 .6 
3460  READ  K(«> 

8470  T*<  Tsteam+273 . 15J/647 .3 

9480  Sum*0 

8490  FOR  N-0  TO  4 

8500  Sum-Sum+K<N)«<1-T>'<NH) 

8510  NEXT  N 

3520  Br*Sum/<  T«<1+K<5)»<1-T)+K<6)*<1-T)  2) )- < 1  - T> 7<K< 7 >*< 1 -T)  2+K<8>> 

8530  Pr-EXP<Br) 

3540  P-22120000-Pr 

8550  RETURN  P 
8560  FNEND 
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APPENDIX  F 

AN  EXAMPLE  OF  REPRESENTATIVE  DATA  RUN 


Month,  date  and  time  :  1 2 : 03 :  1 6  : 5.0  : 24 

NOTE:  Program  name  :  DRP 
Disk  number  *  02 
Old  file  name:  NH05 

This  data  set  taken  on  :  11:02:10:36:28 
Thermocouples  were  defective  at  locations:  0  0 

Tube  Type  is  :  2 

I’ata  Set  Number  -  I  Bulk  Oil  -  0.0  11:02:11:01:06 

T  C  No :  1  2  3  a  5  g  7  o 

Temo  :  70.05  55.07  61.00  70.23  76\68  60.2?  76.21  53.20 

Tina  Tlxgd  Thetab  Tvapr  Peat  Tsump 

56.34  10.03  5.S60E+01  10.01  3.86  -13.0 

Htube  Odp 

».545E*03  8 . 588E  *04 

Data  Set  Number  «  2  Bulk  Oil  7.  -  0.0  11:02:11:30:52 

TC  No:  1  2  3  4  5  6  7  fi 

Temp  :  52.64  42.70  45.49  49.35  58.26  44.85  57.78  47.04 

Twa  Tliad  Thetab  Tvapr  Psat  Tsump 

49.84  10.07  3.939E*01  10.30  3.88  -12.0 

Htube  Qdp 

I  .  I52E*03  4 . S38E*04 


Data  Set  Number  •  3  Bulk  Oil  7.  -  0.0  11:02:11:38:07 


TC  No:  1  2  3  4  5 

Temp  :  34.63  30.25  31.40  33.49  37*59 

Iwa  Til qd  Thetab  Ivapr  Psat 

33.72  10.00  2 . 357E+0 1  10.69  3.83 

Htube  Odp 

7 . 420E  *02  1.749E*04 


6 

31.38 

Tsump 

-11.9 


38.00 


8 

32 . 96 


Data  Set  Humber  -  4  Bulk  Oil  7.  -  0.0  11:02:1  1:49:56 


TC  No:  1  ?  3  4  6 

Temp  :  27.20  24.92  25.40  26.60  29.04 

Twa  Tliad  Thetab  Tvapr  Psat 

26.78  10.07  1.663E*01  11.17  3.89 

Htube  Qdp 

5.2ME  +  02  8 . 664E*03 


6 

25 . 25 
T  sump 

-11.9 


7 

29.30 


8 

26.40 


Data  Set  Number 


5  Bulk  Oil  X 


11:02:1! : 57 : 02 


T  C  No :  1  2  3  4  5  6  7  8 

Temp  :  22.10  21.03  21.06  2' .41  2?! M  21.55  23.11  21.65 

Twa  Tlnd  Thetab  Tvapr  Psat  Tsump 

21.89  9.95  l.189E*01  I '.47  .3.90  -12.1 

Htube  Odp 

3.0946*02  3.680E*03 


Data  Set  Number  *  6  Bulk  0 i  1  '/. 


11  :02:  12:05:46 
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imwi.m-Ui u nmm  i-tu  ueuj nu  «.«  *u  ii  m.i> 


TC  No: 
Temp  : 
Tuia 
17.95 
Htube 


1 

18,14 
T1  iqd 
10.19 
Qdp 


2  8 
17.58  17.82 

The  tab 
7.737E+00 


. 728E+02  ! . 337E  *09 


6  5  8 

17.57  19. IE  18.03 

Tvapr  Psat  Tsump 
12.00  3.96  -12.2 


7  8 

18.43  16.84 


Data  Set  Number  *  7  Bulk  Oil  7.  -  0.0  I  I  :  02  :  1 2 :  1 4  :  1 8 

TC  No:  1  2  3  4  5  6  7  8 

Temp  :  15.58  15.08  15.14  14.87  15. 6E  15,04  15.46  14.55 

Twa  Tliqd  Thetab  Tvapr  Psat  Tsump 

15.17  10.12  5 . 042E+00  12.34  3.92  -12.2 

Htube  Qdp 

1 .221E+02  6 .  I56E+02 

NOTE:  07  X-Y  pairs  were  stored  m  plot  data  file  P05 
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APPENDIX  G 

LISTING  OF  CALIBRATION  COMPUTER  PROGRAM  (TCAL) 


trio  ♦  F  It  E  NAME:  TCAL 

MO  !  REVISED:  October  4.  1984 

ICO  » 

I  CO  COM  /Cc/  C<  7 ) 

140  DIH  En f ( 4  ) .  T ( 4 ) . D< 4 > 

ISO  DATA  0. 1008609! .25727.90369.-787345.8295.78029995.31 

160  DATA  -9247486589. 6. 97688E11 .-2. 66192E13. 3. 94078E14 

1  70  READ  C<*> 

130  CLEAR  709 

190  BEEP 

200  INPUT  "ENTER  MONTH.  DATE  AND  TIME  < MM : DD : HH : MM : SS) " . B5 
210  J*P 

220  OUTPUT  709  :"TD" :BS 

230  OUTPUT  709 ; "TD” 

24 0  ENTER  709 : AS 

241  PEEP 

242  INPUT  "WANT  A  HARD  COPY  ( 1 *Y  ,0-N)7" . lhp 

244  IF  Ihp-1  THEN  PRINTER  IS  701 

245  F  Ihp-0  THEN  PRINTER  IS  1 

250  PRINT  USING  "1  OX , ""Month .  date  and  tine  -  "”.14A":A$ 
260  BEEP 

270  INPUT  "ENTER  INPUT  MODE  (1-3054A.  2-FILE>". Im 

230  IF  In-1  THEN 

290  BEEP 

300  INPUT  "GIVE  A  NAME  FOR  DATA  FILE",D  file? 

310  CREATE  BDAT  D_flleS,20 

320  ELSE 

330  BEEP 

340  INPUT  "GIVE  NAME  OF  EXISTING  F ILE" . D_f I  I e$ 

350  BEEP 

360  INPUT  "ENTER  NUMBER  OF  DATA  RUNS  STORED". Nrun 

361  BEEP 

362  INPUT  "GIVE  A  NAME  FOR  OUTPUT  FILE",OfiIeS 

363  CREATE  BDAT  Ofil-S.IO 

364  ASSIGN  «TiIeo  TO  OflleS 

365  K* 1 

370  END  IF 

380  ASSIGN  9F i 1 e  TO  D_TiIe$ 

390  IF  Im* 1  THEN 

400  BEEP 

410  *  INPUT  "ENTER  BATH  TEMPERATURE" . T_bath 

411  OUTPUT  7 1 3 : ” T 2R2E" 

412  WAIT  2 

413  ENTER  71 3 :  T  bath 

420  OUTPUT  709 : "AR  AF29  AL33  VRI" 

430  EOR  1-0  TO  4 

440  OUTPUT  709: "AS  SA" 

450  ENTER  709 : EnT  < I ) 

460  NEXT  I 

470  ELSE 

480  ENTER  ■Eile:T_bath,Eraf<*> 

490  END  IE 

500  J-JH 

510  FOR  1-0  TO  4 

570  T  <  I  )-FNT-/sv<  ABStEnf  ( I )  >  ) 

530  D< I > - T  bath- T <  I  > 

540  NEXT  I' 

541  IF  K-1  THEN 

542  FOR  L  *0  TO  4 

544  OUTPUT  «FHeo:T<L).D(U 
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!'■  U‘»  J!«  l»*J 


T*  ms1 


505 
50(5 
50y 
550 
560 
570 
530 
530 
600 
610 
620 
630 
600 
660 
660 
670 
630 
630 
700 
710 
720 
730 
7  00 


NEX  1  L 
Gnrn  ?oo 

END  IF 
PRINT 

PR  III  r  US  IMG 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
IF  Im-2  THEN  700 
BEEP 

INPUT  "OK  ro  ACCEPT  THIS  SET  < 1 -Y , 0»N) T* ,0k s 

IF  Oks-O  THEN 

J-J-1 


'  1  OX  ,  ■■■' th  temperature  *  '*** .  3D  -  3D . ' 

1 1  OX  Thermocouple  readinqs  (  Deq  C>: . 

'  MX .  5  ( IP . DO  ,  3X  )  .  ,?5X"  ;  T  <  »  ) 

’  1  OX . "“D isc repanc i es  (Deq  C): . 

M  I  X , 5  <  M2 . DO . OX ) , 24X“ ; D( • ) 


Gnrn  300 

ELSE 


OUTPUT  “*File:T  bath ,£nf < • ) 

END  IF 
DEEP 

IF  I  m  *  1  T  HF  M 

51IPUI  "14 ILL  THERE  BE  ANOTHER  SET  (1  *Y  ,0 .Go  on 

IF  Go_on- 1  THEM  330 

ELSE 


(Deq  C ) 


; T  _bath 


7S0  IF  J<Nrun  THEN  330 

760  END  IF 

770  PRINT 

780  IF  Im*1  THEM 

730  PRINT  USING  "I  OX."' 

J  ileS 

300  ELSE 

310  PRINT  USING  "MX."' 
fileS 

3M  PRINT  USING  "MX."' 

820  END  IF 

320  ASSIGN  *File  TO  * 

808  END 

360  (IFF  F NTvs>/< Emf  > 

860  COM  /Cc/  C <  7  » 

3  70  $um»C<0) 

880  FOR  Ml  TO  7 
830  Sum-Sum »C< I >*Emf  I 
9H0  NEXT  I 
3M  RETURN  Sum 
320  FNEND 


NOTE:  "".OD,""  data  sets  are  stored  in  file  "".K|A";J,D 

NOTE:  Above  analysis  was  performed  from  file  ’"',MA":D_ 
Output  data  are  stored  in  Tile  ""  .  IW'  ;0f  i  le? 
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APPENDIX  H 
SAMPLE  CALCULATION 


Data  run  number  5  (saturation  temperature  was  10  °C  and 
heat  flux  was  about  20  kW/m2)  was  chosen  for  the  sample 
calculation. 

A.  TEST- SECTION  DIMENSIONS 


Do 

= 

0.01588 

(m) 

Di 

= 

0.01270 

(m) 

Di 

= 

0.01143 

(m) 

= 

0.01588 

(m) 

L 

= 

0.20320 

(m) 

Lu 

0.07620 

(m) 

MEASURED 

PARAMETERS 

vs 

= 

3.58 

(volts ) 

Is 

= 

2.27 

(volts ) 

T  l 

= 

34.69 

(°c) 

t2 

= 

30.25 

(°C) 

T3 

= 

31.40 

(°C) 

t4 

= 

33.49 

(°C) 

T  5 

37.59 

(°C) 

t6 

= 

31.38 

(°C) 

t7 

r 

38.00 

(°C) 

Ta 

= 

32.96 

(°C) 

I  s  at 

= 

10.00 

(°C) 

k  c 

= 

344.00 

(W/m.K) 
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C.  OUTER  WALL  TEMPERATURE  OF  THE  BOILING  TUBE 


p  =  ir  £  =  3.1416  x  0.01588  =  0.0499  (m) 

Ac  =  71  (D20  -  D*i)/4  =  3.1416  x  (0.015882  -  0.01272)/4 
A  c  =  714  x  10  "7  (m2 ) 

V  =  25VS  =  25  x  3.58  =  89.50  (volts) 

I  =  Ig  =  2.27  (amp) 

Q  =  VI  =  89.50  x  2.27  =  203.17  (W) 

8 

Tav<f  (  £  )  /  8  *  33.72  (°C) 

n  =  1 

\o-  UV  ln<V  Dl>/2  *L  kc> 

=  33.72  -  203. 17(ln(15. 88/11. 43)/2x3. 1416x203. 2x344) 

\0  =  33.57  ( °C ) 

'  ^wo 

0  =  Tv,-  T  =  33.57  -  10.00  =  23.57  (°C) 

D  sat  ' 

D.  PROPERTIES  OF  R-114  AT  FILM  TEMPERATURE  (REF.  24},  {REF. 
25} 


T  ,  =  (f  +  T  J/2  =  21.78  (°C)  =  294.78  (K) 

f  v  wo  sat' '  v 

log  U  (10  -JNs/m2  )  =  -4.4636  +  1011.47/Tf 
M  =  exp(-4.4636  +  (1011.47/294.78))  x  10  '3 
U  =  356  x  10  ■*  (Ns/m2) 

P  =  36.32  +  61.14  j1'3  +  16.42  j  +  17.48  j1/2  + 

1.12  j2  (lb/ft3) 
j  =  1  -  T,  /T 

Tc  =  Critical  Temperature  =  753.95  (R) 

T  =  294.78  (K)  =  530.6  (R) 

j  =  1  -  530.60/753.95  =  0.296 
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P  =  91.55  (lb/ft3)  =  1466.49  (kg/m3) 

v  =  u  /  P  =  2.428  x  10  -7  (m*/s) 

k  (W/m.K)  =  0.0710  -  0.000261  Tf  (T^  in  °C) 

k  =  0.0653  (W/m.K) 

Cp  (kJ/kg.K)  =  0.4  ♦  1.65  x  10"3Tf  +  1.51  x  10  -*Tf  - 
6.68  x  10  -lflTf3  (Tf  in  K) 

=  1002  ( J/kg.K) 

a  =  k/  Pq>=  4.444  x  10  *■  (m2/s) 

8  =  -(Ip  /  AT)/  P 

8  =  -(1465.95  -  1466.49)  /  ((294.98  -  294.78)  x  1466) 

8  =  0.00184  ( 1/K) 

pr  =  v  /  a  =  5.464 

E.  HEAT -FLUX  CALCULATION 


Average  natural-convection  heat- transfer  coefficient 
at  non-boiling  ends: 


hp  '1/2Luf  1/6 


\ 


/  0.60  +  0.387 


II  +  (0.559/Pr)9/16)8/27 


h  =  282  (W/m’.K) 


Heat- transfer  rate  through  non-boiling  ends: 


Q  =  (h  p  kcAc)1/I  0bTanh(mLu) 
m  =  (hp/kcAc)l/J  =  23.94  (1/m) 

=  13.15  (W) 

QLos=  2  x  Qf=  2  x  13.15  =  26.30  (W) 

Heat  flux  through  active  boiling  surface: 

Q  =  Qh'  ^losI  203-17  ‘  26-30  =  176.87  (W) 

q  =  Q  /  A  b 

Ab  =  DqL  =  10.137  x  10  °  (ml ) 
q  =  176.87  /  0.010137  =  17447.96  (W/m2 ) 

h  =  q  /  Gh  =  740.26  (W/m2 . K ) 

The  following  are  the  results  obtained  from  the  computer 
by  running  the  data  reduction  program  (See  Appendix  F). 

q  =  17490  (W/m2) 

^  =  23.57  (°C) 

h  =  742  (W/m2 .K) 


APPENDIX  I 

UNCERTAINTY  ANALYSIS 

The  same  data  set  (run  number  5)  that  was  used  for  the 
sample  calculation  was  chosen  for  the  uncertainty  analysis; 
therefore,  the  dimensions  of  the  test  section,  and  the  meas¬ 
ured  and  calculated  parameters  found  in  the  sample  calcula¬ 
tion  were  used  in  this  analysis.  All  uncertainties  are 
presented  as  a  percentage  of  the  calculated  parameter. 

A.  UNCERTAINTY  IN  SOURCE  HEAT- TRANSFER  RATE 
Qh  =  VI  (W) 

I  =  Ig  =  2.27  (amp)  <$I  =  ±0.025  (amp) 

V  =  25V  ;  V  =  3.58  (volts)  5VS  =  ±0.05  (volts) 

where : 

6  =  uncertainty  in  measurement  and  calculation 

6Qh/Qh  =  ((5vs/vs)2  +  (Sl/I)2)1  /2 

5Qh/Qh  =  ((0.05/3. 58 )z  +  (0 . 025/2 . 27 ) 2 ) 1  '  2 

=  1.80  percent 

B.  UNCERTAINTY  IN  SURFACE  AREA 

Ab  =  ^oL 

DQ  =  15.88  (mm)  6DQ  =  0.1  (mm) 

L  =  203.20  (mm)  $L  =  0.1  (mm) 

^/A^  =  ((6D0/D0)2  +  (  6L/L) 2  ) 1  '  2 

6Ab/Ab  =  ( (0 . 1/ 15 . 88 ) 2  +  (0. 1/203. 2)2)1  12 

=  0.63  percent 


C.  UNCERTAINTY  IN  WALL  SUPERHEAT 


w. 

k 


m 

'A 


i 


Lv,' 

Krv 

lV 


:• 


AT  = 


sat 


T  =  10.00  (°C) 

Sdt 


5T  =  0.5  (°C) 

X,o  =  ^vg-  (ln(r>2  /Dx  )  /  2  tt  L  kc)  . 


r  = ( ) 

avg  \  4_i,  n' 
n  =  l 


/  8 


where : 

TJj  =  thermocouple  readings  (See  Appendix  G) 


W  33.72  (°C) 


S.D.  = 


wi)/7  ),m 

n=»  1 


2.687  (°C) 


where : 

S.D.  =  standard  deviation 
Since  logarithmic  term  in  equation  of 
when  compared  to  standard  deviation,  this  term  can  be 
neglected  for  uncertainty  analysis  (i.e.,  \Q=  Tsat). 


"Iwo"  is  too  small 


f  =  33.72  (°C) 

wo 


6T 


wO 


2  x  S.D.  =  5.37  (°C) 


AT  = 

T  -  T  =  33.72 

wo  sat 

<5A  T/A  T 

=  ((6  WAT)  2  +  ( -6  ' 

6 A  T/AT 

=  ((5. 37/23.72)*  «■ 

<SA  T/AT 

=  22.7  percent 

D.  UNCERTAINTY  IN  HEAT  FLUX 

q  = 

<V  2Qp)/*b 

%  - 

203.17  (W) 

Assuming 

the  same  proportion  : 

%  - 

13.15  (W) 

%  - 

2Qf  =  176.87  (W) 

6Q 


3.65  (W) 


SQ  =  0.23  (W) 

F 
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6q/q  =  (  (<SQ/Qh-  2QfF  ♦  2(^1*  - 

(64/Afc)2  )2'2 

6q/q  =  (  (3. 65/176. 87)2  ♦  (0 . 46/ 176 . 87 ) 2  + 

( -0 . 0063 ) 2  )2 '2 
=  2.17  percent 

UNCERTANITY  IN  BOILING  HEAT-TRANSFER  COEFFICIENT 


h  =  q  /  A  T 

6h  /  h  =  (  ( ($q  /  q ) 2  +  ( -6  A  T/A  T)  2  )‘  '2 

5h  /  h  =  (  (0 . 0217 ) 2  +  ( -0 . 227  )  2  j1'2 
=  22.8  percent 
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